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SUMMARY 


This  report  contains  a  comprehensive  description  of  a  computer  program  developed  for 
the  numerical  evaluation  of  the  helicopter  noise  equations  derived  in  USAAVLABS 
TR 68-60,  "Studies  of  Helicopter  Rotor  Noise."  It  is  completely  self-contained  in 
that  the  program  details  aie  described,  starting  from  two  basic  acoustic  equations  and 
covering  methods  by  which  these  equations  are  applied  to  the  rotor  noise  problem. 
Program  flow  diagrams  and  a  complete  listing  are  presented  together  with  input  instruc¬ 
tions  and  sample  inputs  and  outputs.  The  program  is  written  in  FORTRAN  IV  for  the 
CDC  3300  Computer. 
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1.0  INTRODUCTION 


The  rotor  sound  field  is  calculated  according  to  the  general  acoustic  equation  for  aero¬ 
dynamic  forces  in  motion  which  is  derived  in  Section  3.3  of  USAAVLABS  TR  68-60: 
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This  equation  is  written  in  tensor  notation,  where  the  i  denotes  summation  over  the 
three  component  directions,  and  the  brackets  denote  evaluation  at  the  appropriate 
retarded  time  t‘  =  t  -  s/c  which  is  the  time  at  which  the  source  generated  the  sound 
reaching  the  observer  at  time  t. 

p(t)  -  p  is  the  instantaneous  acoustic  pressure. 


x,  are  the  coordinates  of  the  observer. 
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y.  are  the  coordinates  of  the  source. 
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M  is  the  component  of  the  source  Mach  number  in  the  direction  of  the 

s  L 

observer . 
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is  the  distance  between  source  and  observer. 


F,  are  the  components  of  the  aerodynamic  force, 

c  is  the  speed  of  sound  . 


In  addition,  the  near  field  pressure  fluctuations  are  calculated  according  to  Low  son  ^ 
as  follows: 
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1  .  Lowson,  M.V.,  "The  Sound  Field  for  Singularities  in  Motion",  Proceedings  of  the 
Royal  Society,  Volume  A286,  pp.  559-572  (1965).  (Equation  18). 
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and 

M. 

is  the  source  Mach  number 
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is  the  component  of  M  in  the  i-direction. 

I 

In  all  cases  the  sum  of  Equations  (1)  and  (2)  is  calculated,  but  as  shown  in 
Section  6  of  USAAVLABS  TR  68-60  the  component  due  to  (2)  is  negligible  at  any 
significant  distance  from  the  rotor. 

Equations  (1)  and  (2)  are  solved  by  a  numerical  method  which  has  been  programmed 
for  digital  computation.  This  solution  is  exact  to  the  extent  that  no  approximations 
are  made.  The  accuracy  of  the  solution  is  only  limited  by  that  with  which  the  real 
loads  and  motions  experienced  by  the  rotor  dynamic  system  can  be  represented  by  the 
model  used.  The  sound  field  for  the  simulated  system  is  calculated  accurately  at  any 
point  in  the  near  or  far  field . 

The  sound  generated  by  a  rotor  blade  in  motion  is  the  result  of  the  distributed  aero¬ 
dynamic  pressure  acting  over  its  entire  surface.  Rotational  noise,  which  is  the  subject 
of  this  study,  is  defined  as  that  component  of  the  sound  field  which  is  directly  attrib¬ 
utable  to  the  lift  and  drag  forces  acting  on  the  blade.  Strictly,  the  entire  spanwise 
and  chordwise  distributions  of  these  components  should  be  taken  into  account,  but  for 
the  sake  of  numerical  expediency  it  is  necessary  to  simulate  the  actual  distributions 
by  a  discrete  set  of  point  loads.  The  implications  of  doing  so  are  fully  discussed  in 
Section  5  of  USAAVLABS  TR  68-60 .  The  spanwise 
loading  distributions  are  divided  into  a 
number  of  segments,  each  of  which  is 
represented  by  two  single  force  components, 
lift  and  drag.  This  model  is  illustrated  in 
Figure  1 .  The  point  of  application  of  each 
force  pair  then  becomes  an  acoustic  source 
which  generates  sound  according  to  Equations 
(1)  and  (2).  The  F-  in  those  equations 

are  the  forces  acting  upon  the  air  and  are 
therefore  opposed  to  the  lift  and  drag  forces 
acting  upon  thu  blade. 

If  an  arbitrary  set  of  orthogonal  axes  X,Y,Z 
are  defined.  Equations  (1)  and  (2)  can 
be  combined.  Using  also  a  more 
convenient  notation, 

Figure  1.  Discrete  Representation 
of  Aerodynamic  Loading. 
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Ap(t)  is  the  observed  sound  pressure  at  time  t  due  to  the  aerodynamic  forces  acting  at 
any  particular  point  on  the  blade  whose  components  in  the  X,Y  and  Z  directions. 


at  the  retarded  time  t  -  — ,  were  Fv,  F.,  and  F.,. 

c  X  Y  Z 


The  dot  denotes  differentiation 


with  respect  to  time.  The  quantities  x,  y  and  z  are  the  coordinates  of  the  observer, 
relative  to  the  source  in  the  X,  Y  and  Z  directions.  The  negative  sign  accounts  for 
the  use  of  the  blade  loads  which  are  equal  and  opposite  to  the  forces  which  act  on  the 


It  can  be  seen  from  this  equation  that  the  force  generates  sound  in  two  ways:  through 
its  fluctuation  in  time  and  its  fluctuations  in  position  which  cause  its  accelerations 
toward  the  observer  (M^) .  The  term  (1  -  M^)  in  the  various  denominators  is  essen¬ 
tially  a  Doppler  effect  which  amplifies  the  sound  radiated  in  the  direction  of  motion. 

Dipole  sound  (the  term  in  F  )  is  amplified  by  the  factor  (1  -  M  )  and  quadrupole 

-3  s 

sound  (the  accelerative  term  in  F  •  by  (1  -  .  Thus  the  second  term  becomes 

increasingly  important  as  the  velocity  in  the  direction  of  the  observer  increases. 

The  elements  of  the  method  of  computing  the  sound  fieid  of  a  complete  rotor  are  as 
follows: 

(1)  Define  *he  geometry,  blade  loading  and  resulting  blade  motions  of  the 
rotor  as  a  function  of  time. 
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(2)  Define  the  position,  attitude  and  velocity  of  the  rotor  with  respect  to 
the  observer  at  the  time  t  . 

(3)  Calculate  the  orientation  and  magnitude  of  each  of  the  elemental  blade 
airloads  at  its  appropriate  retarded  time. 

(4)  Calculate  and  integrate  the  observed  sound  pressures  due  to  all  the  s>  $. 

(5)  Repeat  operations  (2)  through  (4)  for  a  series  of  successive  time  int  s 
t  to  construct  a  time  history  of  the  acoustic  pressure  amplitude  and 
harmonically  analyze  this  to  obtain  its  frequency  spectrum. 

A  computer  program  has  been  written  to  perform  these  operations  taking  account  of  the 
following  variables: 

(a)  Number  of  rotors  . 


(b)  Number  of  blades. 

(o')  Rotor  diameter . 

(d)  Helicopter  position  in  space  with  regard  to  a  fixed  observer 
including: 

1.  Rotor  hub  vertical  displacement. 

2.  Rotor  hub  horizontal  displacement. 

3.  Rotor  shaft  inclination. 

(e)  Rotor  hub  motion  including: 

1 .  Linear  velocity . 

2.  Roll  angular  velocity. 

3.  Pitch  angular  velocity . 

(f)  Rotor  blade  tip  speed . 

(g)  Articulated  blade  motions. 

(h)  Blade  motion  resulting  from  flapping  bending,  edgewise  bending, 
and  twist  about  the  elastic  axis. 
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(i)  Phase  relationships  between  rotor  angular  positions  for  multirotor 
vehicles. 

(j)  Time  dependent  blade  loading  for  at  least  20  radial  blade  segments. 
The  blade  loading  includes  loading  normal  to  the  hub  plane,  radial 
loading  in  the  hub  plane,  and  tangential  loading  in  the  hub  p  une. 
These  loadings  are  treated  as  loads  at  a  chordwise  point  on  the  blade 
element  with  direction  determined  by  orientation  of  the  blade 
element  in  space.  Variation  in  disc  loading  is  considered  as  a 
function  of  variation  in  blade  loading  with  the  number  of  blades 
and  blade  chord  held  constant. 

This  program,  code-named  HERON  1,  is  described  in  detail  in  Sections  4  and  5. 

The  axis  transformations  and  outlines  of  the  computational  steps  are  discussed  in  the 
following  sections. 
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2.0  COORDINATE  SYSTEMS  AND  TRANSFORMATIONS 


2.1  Rotor  Axes  x1  ,y* ,z‘ 

Since  blade  motions  are  measured  relative  to  the  rotor  shaft  it  is  desirable,  for  the 
direct  utilization  of  experimental  data,  to  specify  all  blade  responses,  including 
flapping,  lagging  and  elastic  deformations  with  respect  of  a  set  of  "shaft  axes". 


Figure  2.  Shaft  Axes  ,y' ,z')  and  Flight  Path 
Coordinates  (x,y,z). 

Accordingly,  a  set  of  right-handed  orthogonal  axes  x',  y',  z'  are  chosen  with  the  z' 
axis  coincident  with  the  rotor  shaft  and  positive  downward  (Figure  2) .  The  longitu¬ 
dinal  axis  x'  lies  in  the  vertical  plane  through  the  aircraft  velocity  vector  V  . 

2.2  Flight  Path  Axes  x,y, z 

This  system  is  introduced  merely  to  simplify  the  transformations  between  the  x' ,y' ,z' 
and  X,  Y,  Z  systems.  The  origin  of  this  system  is  coincident  with  that  of  the  rotor 
axes.  However,  the  x  y  plane  is  horizontal  with  the  x  axis  in  the  vertical  V-x1 
plane.  The  aircraft  rotation  is  specified  with  respect  to  this  system  in  the  following 
sequence  (Figure  3). 
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Figure  3.  Pitch  and  Roll  Rotations. 

(1)  The  pitch  angle  0  is  measured  clockwise  about  the  positive  y  axis, 
rotating  x  to  x'. 

(2)  The  roll  angle  is  then  measured  clockwise  about  the  x'  axis,  rotating 
y  to  y'  and  z  to  z'. 

(Vehicle  angular  displacements  are  measured  in  pitch  and  roll  only,since  rotor  yaw  is 
simply  a  shift  in  azimuth  reference.) 

The  shaft  inclinations  to  the  vertical  z  axis,  measured  in  the  x  z  and  y  z  planes, 
are  91  and  <}>'  so  that  (Figure  4)  9  =  0'  and  tan  <}>  =  cos  9  tan  <$>' . 


Figure  4.  Shaft  Inclinations. 
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The  transformations  between  shaft  and  flight  path  axes  are: 

x  =  x‘  cos  9  +  y1  sin  <J>  sin  0  +  z1  cos  <{>  sin  0 

y  =  y1  cos  <J>  -  z '  sin  <|>  (4) 

z  =  x‘  sin  9  +  y'  sin  $  cos  0  +  z'  cos  $  cos  9. 

2.3  Fixed  Axes  X,Y,Z 

This  coordinate  system  is  used  to  define  the  aircraft  and  observer  positions  in  space. 
The  axes  are  fixed  in  space  with  arbitrary  origin,  and  XY  is  the  horizontal  ground 
plane.  Z  is  measured  positive  vertically  upwards. 

The  X  Y  and  x  y  planes  are  thus  parallel,  and  the  angle  between  the  x  and  X 
axes  is  defined  as  X  . 


x 


Figure  5.  Fixed  Axes  (Plan  View). 

Thus,  if  the  rotor  hub  coordinates  are  X^,  Y^,  Z^,  any  point  in  the  X,  Y,  Z 
system  is  defined  by 

X  =  X^  +  x  cos  X  +  y  sin  X 
Y  =  Yj_j  +  x  sin  X  -  y  cos  X 

z  =  ZH  '  2  (5) 


8 


The  coordinates  of  an  observer  at  X0,Yq/Zo  ,  with  respect  to  the  source  at  x,y,z, 
are 


x  =  X  -X 
o 

y  =  Y«-Y 


The  distance  s  between  source  and  observer  is 


s 


2  1 
1  2 


y  +  5  ) 


2.4  Blade  Element  Displacements 


(7) 


The  x1,  y' ,  z‘  coordinates  of  each  specified  blade  loading  point  are  calculated  as  a 
function  of  azimuth  angle,  flapping  and  lagging  angles  (where  applicable),  and 
normal  and  in-plane  elastic  displacements.  The  rotor  azimuth  angle  is  measured  from 
the  negative  x1  axis,  clockwise  about  the  positive  z'  axis,  qj  =  Qt. 


Figure  6.  Blade  Element  Displacements  in 
Rotor  Axes. 
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Referring  to  Figure  6, 


is  the  flap  hinge  offset  from  the  shaft  centerline, 
is  the  lag  hinge  offset . 

C  is  the  lag  angle  projection  on  the  x*  /'  plane  (positive  in  the  direction 
of  rotation)  . 

P  is  the  flap  angle  between  the  undeformed  blade  axis  and  the  x'y'  plane, 
positive  in  the  negative  z'  direction. 


r\j(r)  is  the  elastic  displacement  of  blade  station  r  parallel  to  the  x'y'  plane, 
normal  to  the  undeformed  blade  axis  and  positive  in  the  direction  of 
rotation. 


n 


is  the  elastic  displacement  normal  to  the  undeformed  blade  axis,  in  the 
plane  containing  the  z'  axis  and  positive  in  the  negative  z'  direction. 


These  coordinates  thus  take  account  of  all  possible  motion  of  the  blade  axis  upon  which 
the  point  aerodynamic  loads  are  assumed  to  act. 

If  the  flapping  hinge  is  outboard  of  the  lagging  hinge,  the  coordinates  of  the  radial 
station  r  are: 

x'  =  -  cos  ip  -  j(r  -  e^)  cos  f  +  (e^.  -  e^)  j  cos  (ip  +  C) 

+  q  y  sin  (t|j  +  C)  +  r)  sin  0  cos  (ij>  +  £) 

+  r| y  cos  (i|!  +  £)  +  r|  ^  sin  0  sin  (<|i  +  £) 
z1  =  -  (r  -  ef)  sin  p  -  nN  cos  p  (8a) 

If  the  flapping  hinge  is  inboard  of  the  lagging  hinge,  we  have 


- 


y  =  - 


(£  sin  i|i 


Cf> 


cos 


P  +  («, 
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x  =  - 


y  = 


(e^.  +  cos  p)  cos  ij)  -  (r  -  e^)  cos  p  cos  (iji  +  C) 
+  q  y  sin  (ij)  +  £)  +  q  ^  sin  p  cos  (ijj  +  £) 

(e^  +  cos  p)  sin  iji  -  (r  -  e^)  cos  P  sin  (tj)  +  0 
-  q^.  cos  (iji  +  (>  )  +  q  ^  sin  P  sin  (iji  +  £) 


z'  =  -  (r  -  ef)  sin  p  -  qN  cos  p 


(8b) 


It  is  now  possible,  using  Equations  (4),  (5)  and  (8),  to  find  the  fixed  coordinates 
X,  Y,  Z  of  any  blade  station,  defined  in  convenient  rotor  coordinates  as  a  function 
of  time,  and  subsequently,  using  Equations  (6)  and  (7),  the  displacement  components 
x,  y,  z  and  s  of  the  sound  Equation  (3) . 

2.5  Blade  Element  Velocities  and  Accelerations 

We  also  require  the  velocity  and  acceleration  components  defined  in  the  sound  Equation 
(3)  as  M,  Ms  and  M$  which  are  the  source  Mach  number  and  rate  of  change  of 

Mach  number  in  the  direction  of  the  observer.  Resolving  the  source  Mach  number  into 
its  three  components  we  have 


M 

II 

+ 

H*. 

~r 

Mz ' 

(9) 

M 

s 

=  —  M  +  1-  M 
s  X  s  Y 

+  -  M7 
s  Z 

(10) 

and 

M 

s 

-  1  M  +  1  M 
s  X  s  Y 

+  -  m7 

s  Z 

(11) 

where  My,Mv,  and  My  are  simply  — ,  —  and  — •  respectively;  that  is,  the 
AY  Z.  c  c  c 

velocity  components  of  the  source,  relative  to  the  stationary  air,  expressed  as  a 

fraction  of  the  atmospheric  speed  of  sound.  To  obtain  the  first  and  second  time 

derivatives  of  X,Y,  and  Z  it  is  necessary  to  differentiate  Equations  (5),  and 

subsequently  Equations  (4)  and  (8).  From  Equation  (5), 
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+  *  cos  x  "  x  X  s'n  X  +  y  sin  x  +  Y  X  cos  X 


X 

Y 


Z 


*  =  *H 

=  X|_|  +  (x  +  y  x)  cos  x  +  (y  -  X  x)  sin  x 

•  • 

Y  =  Y|_|  +  x  sin  x  +  x  X  cos  X  "  y  cos  X  +  y  x  sin  x 

=  Y|_j  -  (y  -  X  x)  cos  x  +  (x  +  y  id  sin  X 

z  =  ZH  -  z  (12) 

XH  +  (x  +  2y  X  -  X  X2  +  y  x)  cos  x  +  (y  -  2x  x  -  yx2-  xj()  sin  X 

Yh  -  (y  -  2  xx  -  yx2  -  xx)  cos  x  +  (x  +  2yy  -  xy2  +  y  x)  sin  X 

ZH  +  z  (13) 


and  from  Equation  (3), 


x  =  x'  cos  8  -  x'  9  sin  8  -  (z'^>  -  y1)  sin  (j>  sin  8 

+  (z1  +  y'<j>)  cos  <f>  sin  8  +  y'  8  sin  <{>  cos  8  +  z'  8  cos  $  cos  8 

y  =  (y1  -  z'  $)  cos  ({>  -  (y 1  q>  +  z)  sin  <£ 

z  =  z  -  x'  sin  0  -  x*  0  cos  8  +  (y1  -  z '<^>)  sin  q>  cos  6 

+  (y'<j>  +  z1)  cos  cos  8  -  y'0  sin  q>  sin  0  —  z'8  cos  <{>  sin  0  (14) 


x  =  (x1  -  x’92)  cos  0  -  (2x'0  +  x'  8)  sin  0 
+  (2y 1  <j>0  +  2  z'0  +  z'8)  cos  q>  cos  9 
+  (y1  -  y'  +  02)  -  2z'ci>  -z'q>)  sin  q>  sin  0 
+  (2y'<j>  +  y'<£  +  z1  -  z1  (<j>2  +  82))cos  <{>  sin  8 
+  (2y'  8  +  y'0  -  2z'<j>9)  sin  $  cos  8 
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y  =  (y‘  +  y’4>2  -  z*$)  cos  $  -  (z‘  -  z'q>2  +  y'$)  sin  $ 

z  =  -  (2x*  0  +  x'  0)  cos  e  -  (x'  -  x'  e2)  sin  0 

+  (2y‘  <j>  +  y1  $  +  z1  -  z1  (cj> 2  +  02)  )  cos  $  sin  0 

-  (2y 1  0  +  y‘  0  -  2  z  <{>  0)  sin  <f>  si  n  0 


-  (2y '  <j>0  +  2z'  0  +  20)  cos  <}>  sin  0 

+  (y1  -  y '  (<j> 2  +  02)  -  2z‘  <j>  -  z1  <|>  -  z'$)  sin  <|>  cos  0  (15) 

Expressions  for  the  derivatives  of  x',  y'  and  z'  are  not  derived  explicitly  for 
reasons  which  will  be  explained  in  Section  3.0. 


2 . 6  Aerodynamic  Forces  and  Their  Derivatives 


The  aerodynamic  force  components  in  the  sound  Equation  (3)  are  derived  from  the 
components  and  F^  which  act  directly  on  the  blade.  F^  is  defined  as  the 

component  of  the  aerodynamic  load  on  the  blade  which  acts  normal  to  the  blade  axis 
(in  its  deformed  condition)  in  the  plane  containing  the  shaft  axis  z' .  It  acts  in  the 
same  sense  as  the  lift  on  the  blade.  F^.  is  the  tangential  component,  again  acting 

normal  to  the  blade  axis  but  parallel  to  the  x'  y'  plane.  Its  sense  is  opposite  to  the 
blade  drag  force.  Referring  to  Figure  6,  we  see  that  the  blade  force  components  in 
the  x‘,  y'  and  z'  directions  are: 


F  ,  =  Fy  sin  (i|i  +  C  )  +  F^  sin  (31  cos  (t|t  +  ? ) 

I  \ 

F  ,  =  -  Fj  cos  (ij>  +  s  )  +  F^|  sin  p1  sin  (tj>+  C  ) 

Fz'  =  ■  FNC0S  P'  (16) 

t 

where  p1  and  C  are  the  blade  slopes  relative  to  the  x'  y'  plane  and  the  radius 
at  azimuth  tjj  respectively.  That  is, 
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The  transformations  which  are  required  to  obtain  the  components  F^,  Fy  and  F^  in 

the  fixed  coordinate  system  are  identical  to  those  used  to  convert  x',  y',  z'  to 

•  *  * 

X,  Y,  Z  (Equations  (4)  and  (5)) .  The  derivatives  F^/  Fy  and  are  also  derived 


from  the  rotor  axis  values  using  similar  transformations  (Equations  (12)  through  (15)), 

•  •  • 

and  again  the  derivatives  in  the  rotor  coordinate  system  F  , ,  F  ,  and  F  •  are  not 

.  x  y  z 

derived  explicitly  for  reasons  outlined  in  the  following  section. 
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3.0  METHOD  OF  SOLUTION  BY  DIGITAL  COMPUTER 


The  technique  for  calculating  the  helicopter  rotor  sound  field  is  best  explained  by  a 
description  of  the  computational  steps  of  the  program  HERON  1 .  Figure  7  is  an 
illustration  of  these  basic  steps  in  flow  diagram  form.  The  program  is  divided  into 
two  phases.  The  *irst  processes  the  input  data  into  a  convenient  form  for  storage  and 
the  second  compiles  the  sound  pressure  level  at  a.  preselected  number  of  field  points. 

3 . 1  Phase  1  Initial  Data  Processing 

For  each  rotor,  the  blade  loading  and  motion  data  is  read  by  the  program  either  as  a 

number  of  arrays,  listing  their  values  at  a  discrete  number  of  points  over  the  rotor 

disc,  or  as  a  series  of  Fourier  coefficients,  one  set  for  each  radial  station.  To  reduce 

the  volume  of  later  computations,  this  initial  phase  of  the  program  calculates  and 

stores  the  displacements  of  each  loading  point  with  respect  to  the  rotor  axes  x',  y', 

z'  according  to  Equation  (7).  The  three  force  components  F  .,  F  .  and  F  ,  are  also 

x  y  z 

using  Equation  (16) .  These  calculations  are  performed,  upon  the  input  data  in  which¬ 
ever  form  it  comes,  at  a  predetermined  number  of  azimuth  stations.  The  history  of 
each  of  these  six  variables  around  the  azimuth  is  then  harmonically  analyzed  for  every 
radial  station,  and  the  Fourier  coefficients  are  stored.  This  operation  is  carried  out  for 
one  blade  only  since  it  is  assumed  that  all  blades  are  loaded  and  respond  in  identical 
manners.  From  the  stored  coefficients  it  is  possible  for  later  elements  of  the  program 
to  interpolate  for  all  required  blade  loading  and  motion  variables  at  any  arbitrary 
azimuth  station . 

3.2  Phase  II  Calculation  of  the  Sound  Field 


Phase  II  of  the  program  is  repeated  in  its  entirety  for  each  required  observer  position, 


3.2.1  "Observer  Time"  t 

The  sound  field  observed  at  the  position  XQ,  Y  ,  Zq  is  calculated  as  a  time  history 

of  the  acoustic  pressure  Ap  at  the  "observer  time"  intervals  0,  At,  2  At,  . . .  T, 
where  T  is  the  fundamental  period.  The  sound  harmonic  amplitudes  are  then  obtained 
through  a  Fourier  analysis  of  this  time  history. 

The  fundamental  period  T,  to  an  observer  moving  with  the  aircraft,  is  the  blade 
passage  period  2tt/  B.  However,  for  a  stationary  observer  this  shifts  by  the  factor 
(1  -  Mq)_1  where  M  's  ^e  oircraft  Mach  number  component  in  the  direction  of  the 

observer. 
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One  approach  is  to  calculate  r  9  observed  period  which  can  then  be  divided  into  an 
appropriate  number  of  intervals.  However,  because  of  the  aircraft  motion,  a  slight 
error  is  introduced  by  this  procedure  which  results  in  the  computation  of  something 
slightly  less  than  a  complete  period.  This  error  increases  with  sound  harmonic  number 
(which  is  equivalent  to  a  reduction  of  wavelength). 

The  alternative  which  has  been  adopted  in  the  program  is  to  calculate  the  time  history 
at  a  point  which  "moves"  with  the  same  velocity  as  the  rotor  hub  but  which  has  the 
average  position  X0,YQ,Z0  during  the  period  T.  The  only  modification  required  for 

the  final  result  is  to  correct  the  observed  frequency  for  the  Doppler  shift  factor 

(1  -  Mq)  .  The  quotation  marks  are  used  because  at  each  time  increment  the  sound 
pressure  is  still  calculated  at  a  stationary  point  which  has  simply  moved  the  same 
distance  as  the  rotor  since  the  previous  time  increment.  All  relative  velocity  effects 
in  the  acoustic  calculation  are  retained. 


For  each  observer  time  t,  the  sound  pressure  Ap  is  calculated  at  the  effective 
observer  positions: 


X1 

0 


X  + 
o 


V'-T> 


Y' 

o 


Yo  +  V  ' 


Z' 

o 


Zo  + 


zh  (t  - 


(18) 


To  calculate  the  Mach  number  component 
Mq  ,  and  hence  the  Doppler  shift  correc¬ 
tion,  it  is  first  necessary  to  calculate  the 
sound  propagation  time  tq  as  follows: 

The  position  of  the  rotor  hub  at  time  t 

is  X  ,  Y  ,  Z  and  the  distance 
H  H  H 

between  the  hub  and  the  observer  at 
Wo-  according  to  Equation  (6), 

is 


Figure  8.  Velocity  Diagram  for  Sound 

Radiation  by  Moving  Aircraft. 
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_  ..2  -2  -2  2 
s  =  lx0  +  yo  +  zo> 


However,  the  sound  arriving  at  the  observer  at  time  t  was  generated  at  the  retarded 
time  t  -  tq  when  the  hub  was  at  X^,Yj_|,Z^  .  The  corresponding  velocity  triangle 

is  shown  in  Figure  8.  The  resultant  aircraft  velocity  is 


V 


+  Y‘ 


IA 


H 


(19) 


and  the  propagation  time  tq  is  given  by 


T0 


(*o*H  +  y0YH  +  £oiH)  +  yoS'H+zo^H)Z  +  s2(c2  '  y2)  ’ 


c2  -  V2 


(20) 

From  the  propagation  time  the  ,  coordinates  at  the  instant  of  sound  emission  can  be 
calculated.  The  component  hub  to  observer  separations  are  Xg,y'Q  and  where 


x1  =  x  +  X. .  T 

0  0  Ho 


=  y'o  +  YH  t0 


zL  =  z  +  Z 


Hto 


(21) 


Hence,  the  required  Mach  number  component  is 


M0  = 


«o*H 

C2  Tn 


(22) 


3.2.2  Retarded  Time 

The  retarded  time  of  each  blade  loading  point  Kjs  to  be  calculated  independently  for 
each  observer  time  instant.  This  requires  solution  of  the  transcendental  equation 
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1 


(23) 


/  2  -2  2 
CT  =  \x  (t  -  t)  +  y  (t  -  t)  +  £  (t  -  t) 

This  equation  can  only  be  solved  by  iteration,  starting  with  some  realistically  chosen 
value  of  t  and  successively  substituting  newly  calculated  values  into  the  right-hand 
side  until  the  solution  converges.  The  right-hand  side  of  Equation  (23)  is  in  fact  a 
very  lengthy  expression  which  includes  Equations  (4),  (5),  (6)  and  (8).  The  following 
procedure  is  followed  to  minimize  the  computational  time  involved. 

For  the  most  inboard  blade  loading  point,  the  rotor  hub  retarded  time  (t  -  tq)  is  used 

as  a  starting  value.  The  RHS  of  Equation  (23)  is  calculated  using  the  first  harmonic 
Fourier  coefficients  of  x',  y1  and  z'  .  As  an  example. 


qr  =  tp  (t  -  t) 

x'  “  a  ,  (r)  -t  a  .  (r)  cos  d>  +  b  .  (r)  sin  tb 

x  o  x,  X1 

The  values  of  y'  and  z'  are  calculated  similarly  and  the  necessary  transformations 
applied  to  compute  x  ,y  and  z.  These  are  substituted  into  Equation  (23)  to 
yield  a  second  approximation  to  the  retarded  time.  The  iteration  proceeds  but  a 
successively  greater  number  of  harmonics  is  admitted  in  the  calculation  of  the  x',  y',  z' 
coordinates  in  successive  iterations.  By  the  time  convergence  is  reached,  the  pre¬ 
determined  limiting  number  of  harmonics  is  admitted.  It  is  found  that  this  method 
considerably  reduces  rhe  volume  of  computation  (through  avoiding  a  large  number  of 
Fourier  summations)  without  increasing  the  number  of  iterative  cycles  to  convergence. 

The  calculations  for  the  next  radial  station  use  the  final  value  of  t  from  the  first 
radial  station,  and  so  on.  Experience  shows  that  an  average  of  four  or  five  iterations 
is  required  for  a  convergent  solution  with  an  error  of  less  than  10  5  seconds. 

3.2.3  Blade  Loads  and  Motions 

The  Fouriei  repiesentation  of  the  load  and  motion  distributions  enables  the  required 
displacements,  velocities,  accelerations,  forces  and  rate  of  change  of  forces  to  be 
computed  with  a  minimum  amount  of  effort.  From  the  retarded  time  t  -  t,  the  rele¬ 
vant  blade  azimuth  angle  q>  is  obtained  and  the  motion  and  load  variables  in  the 
x‘,  y‘,  z1  system  are  syrthesized  as  shown  in  the  following  example: 

i|>  <(>  (t  -  t) 
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x1  ~  axi  ^  axi  cos  k(ji  +  b^,  sinkiji 

k  k 


k  =1 


x'  =  -  O 


K  j  j 

Zka,  sin  k  «j»  -  b  ,  cos  bi|i  . 
I  X  k  X  1 


k  =1 


1 


x'  = 


K  (  j 

N  k2  a  ,  cos  ki|)+b,  sin  k  tj)  • 

f,  l  *k  *k  I 


The  transformations  listed  in  Section  2.  1 .2  are  then  applied  to  derive  all  the 
necessary  components  of  the  sound  Equation  (3)  . 

3.2.4  The  Sound  Field 


The  observed  sound  pressure  increment  generated  by  each  loading  point  is  calculated 
according  to  Equations  (1)  and  (2)  ,  The  inciements  due  to  all  loading  points  on 
all  blades  are  then  summed  ro  give  the  total  pressure  increment  at  time  t.  The  entire 
process  is  repeated  for  successive  time  instants  until  a  history  is  obtained  for  one  com¬ 
plete  (blade  passage)  period.  The  final  step  is  to  Fourier  analyze  this  time  history 

into  its  sine  and  cosine  component  harmonics  (a  and  b  ).  Each  harmonic  sound  pres- 

n  n 

sure  level  is  then  expressed  in  decibels  by  performing  the  transformation 


SPL 


n 


f 


124.6 


) 


dS  re  .0002  p  Bar 


3.2.5  Mu  I  ti  pie  Rotors 


In  its  present  form  the  program  computes  the  sound  harmonics  at  each  field  point  for 
each  rotor  and  outputs  the  results  independently.  If  the  rotors  have  the  same  funda¬ 
mental  blade  passage  frequency  then  the  sound  harmonicsdue  to  the  individual  rotors 
reinforce  each  other  and  the  in-phase  and  out-of-phase  components  are  simply  added 
together  separately.  If  the  rotors  have  different  fundamental  frequencies  (for  example 
the  main  and  tail  rotors  of  a  "single  rotor"  helicopter),  then  the  individual  harmonics 
do  not  interfere  and  the  sound  spectra  are  merely  superimposed.  The  final  block  in 
Figure  7  is  drawn  with  a  broken  line  since  the  program  does  not  perform  this  function 
automatically.  However,  sufficient  information  is  output  (namely  the  frequency, 
cosine  and  sine  harmonic  amplitude  in  lb/ ft2  and  sound  pressure  level  of  each  har¬ 
monic)  to  enable  the  necessary  summations  to  be  performed  rapidly  by  hand. 


22 


4.0  PROGRAM  FLOW  CHARTS,  EQUATIONS  AND  LISTING 


The  computer  program  HERON  1  is  written  in  FORTRAN  IV  for  the  Control  Data  Corpora 
tion's  CDC  3300  utilizing  the  SCOPE  opera*-  .  -ystem .  The  minimum  hardware 
configuration  required  is  as  follows: 

1  .  3300  series  computer  with  32,000-word  core  storage. 

2 .  Card  reader . 

3 .  Line  printer. 

4.  One  scratch  tape  unit  (LUN2) 

A  complete  description  of  the  program  is  contained  in  this  section  which  includes  the 
following  items: 

4.1  A  detailed  flow  chart  for  the  program.  The  noted  equation 
numbers  correspond  to  list  4.2  . 

4.2  A  complete  list  of  programmed  equations. 

4.3  \  program  listing  . 

4.4  A  list  of  program  symbols  and  definitions  (Table  I). 
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List  of  Equations  Computed  by  Program  HERON  1 


Equation  No. 


-  cos  <|>  -  |(r  -  ef)  cos  (3  +  (ef  -  t^)j  cos  (q»  +  C) 
+  r|y  sin  (i|i  +  £)  +  q^  sin  (3  cos  (<{i  +  £) 


b)  y' 


b)  y' 


=  -  sin  tji  -  j(r  -  ef)  cos  p  +  (ef  -  e^)j  sin  (ip  +  C) 

-  r|T  cos  (t|i  +  ?)  +  q  N  sin  (3  sin  (t|>  + 

=  -  (r  -  e  )  sin  P  -  q  N  cos  p  . 

=  -  (e^  +  cos  p)  cos  tf)  -  (r  +  e^)  cos  p  cos  (t|>  +  ^>) 

+  q^.  sin  (i}i  +  C)  +  q  ^  sin  P  cos  (i|)  +  £) 

=  -  (e^.  +  cos  p)  sin  i|i  “  (r  +  cos  P  s'n  (<1>  + 

-  q^  cos  (ijj  +C)  +  q^  sin  p  sin  (t)j  +  C) 


c)  z'  =  -  (r  -  ef)  sin  P  ~  Hn  cos  3- 
3  a)  dqN/dr  =  (qN(n  +  1)  -  qN(n  -  1)^r(n  +  l)-r(n-l)) 

b)  dqN/dr  =  (qN(2)  -  qN(  l))/(r  (2)  -r(l)J 

c)  dq^  /dr  =  (nj(n+  0  "  Hj  (n  “  1))/(r  (n  +  ^  "  r  (n  "  ])) 

d)  dq^-dr  =  ^(2)  -  q^(  l)j/|r  (2)  -  r  ( l)j 

e)  dqN/dr  =  (qN(n)  -  qN(n  -  l))/(r(n)  -r(n-  1)) 

0  dqT/dr  =  (q-j-(n)  -  qy(n  -  l))/(r  (n)  -r(n-l)) 
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4  a)  (31  -  p  +  tan  (dr^/dr) 

b)  C  =  £  +  tan  (dr)^./dr) 

5  a)  F  ,  =  T  sin  (i|j  +  C  )+  sin  P'  cos  (ip  +  C  ) 

b)  F  ,  =  -  F^.  cos  (t|)  +  C )  +  F^j  sin  p'  sin  (q>  +  £  ) 

c)  Fz>  =  -  Fn  cos  P'  • 

M 

6  a)  °on(u)  =  tJtS  “ 

m  =  I 

b)  (u)  =  ^  ^  u  cos  ^(m  -  l)*2trk/M) 

n 

c)  b^  (u)  =  ^  ^  u  sin  |(m  -  1)  •  2irk/Mj 

n 

M 

7  °>  o0n(Fi)  -  *rE  Fi 

m  =  | 

b)  ak  (F  )  =  ^  F;  cos  ((m_  2"k/M) 

n 

c)  bk  =  Ja  S  Fi  S'n  ^  ’  2lxk/M) 

n 

8  a)  xR  =  XH 

b)  =  yh  - Y; 

ZH  =  ZH  Zo 

9  SH  =  V*H  +  +£H 
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10 


V 


12 

a)  x'H 

*H  +  XH  T'« 

b)  y'H  - 

?H  +  Vo 

c)  z-'H  = 

ZH  +  Ti 

13 

M0  - 

(x"HXH  +  9H'i'H  +  5HZH)/a2o 

T1 

0 

14 

=)  X'H  = 

XH  -  v; 

b)  y'h  - 

yh  -  yh  t'o 

c)  - 

ZH  “  T0 

15 

a)  T  = 

2ir/  Q  B 

b)  AT  = 

2h/^M 

16 

a)  Xo  = 

x;w +  xh  ('  ■  t) 

t  =0  (AT)  T 

b)  Y0  = 

Y'„<i'  +  yh('-t) 

Zo  ‘ 

z;(i) +  ZH  (*  -  t) 

17 

♦  o  ’ 

2  IT  I 

=  \  *  T  •  b 

b  =  1  (1)  B 
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a)  41  =  q>ft  +  fit' 
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20 
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b)  xH=  xH  +  xHt* 

0 

c)  VH  -  Yh  ♦  V' 

0 

d>  ZH  =  ZH0  +  V 

e)  0  =  e0  +  0 1* 


f)  4>  =  4>0  +  <j>t' 

g)  <}>  '  =  tan  '  (cos  0  tan  $) 


h)  =  cos2  (<p  cos  0  sec2  <j>  -  0  sin  0  tan  <p) 

a)  w'  =  u'  +  a^  (u‘)  coskifi  +  (u1)  sin  kijj 

b)  u'  =  u'  -  ksi 


c)  u'  = 

a)  x  - 

b)  y 

c)  z 


«•  -  k2  a2 


a^  («')  sin  k  ij>  -fc^  (u1)  cos  k  ijt 
n  n 

(u1)  cos  k  (|i  +  b^  (u1)  sin  k  iji 


x1  cos  0  +  y '  sin  <p‘  sin  0  +  z'  cos  <j>‘  sin  0 
y'  cos  <j >'  -  z'  cos  $>' 

-  x'  sin  0  +  y'  sin  <p'  cos  0  +  z'  cos  <$>'  cos  0 


a)  k  =  x‘  cos  0  -  x1  0  sin  0  -  (z'ty1  -  y')  sin  <|>'  sin  0 

+  (z1  +  y1  <|>')  cos  ($>'  sin  0  +  y'  0  sin  cos  0 
+  z' 0  cos  q>'  cos  0 

b)  y  =  (y‘  -  z‘(j)')  cos  <$>'  -  (y‘<j>'  +  z')  sin  <$>' 
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22 


23 


24 


25 


c)  z 


a)  x 


b)  y  - 


:)  z  =  - 


a) 

b) 

c) 

a) 

b) 

c) 

a) 

b) 

c) 


.  e.  e  cos  6 +  (y' -*'♦') 

+  (yy+i.)  cosV  cose-y'B'si"  V  sin  0-z'8' cosine 
(i,^)cos6-2x’esin0+(2yVe+2i'e)cosV  cos  6^ 

+  (y’-y‘  (i>'2  +  92  )  -  2z,<j’')  sin  t'*'"®  +  (2y  V  +z  -z  '•> 

+  eM)  cost'  sin  e  +  (2y'e  -lz'i'ifi  ™  *' 

+y't'2)  cos  t‘  -(z'  -z'?' 

2x'0  cos  9  -  (x,-*02)  sin  8  +  (2y‘<i>'  +  z' 

_  ^  +02))  cos  t'sine-(2y'e-2z>'0)^^'  »i"8 

-  (2y't'  0  +2z'9)  COS  t'  sin®  +  (y  -y  (i>  +9^ 

-  2  z'4>')  sin  cos  0. 

-  x  -  X. .  -  X  cos  x  -  y  sm  X 
”  0  H 

_  y-Y  -xsinX-y  cos  X 

0  H 

_  7-7  +  z . 

'  Z0  tt 


X  = 

Y  = 

Z  = 

•• 

X  = 

»» 

Y  = 

4» 

z  = 


x  +  X  cos  X  +  y  sin  x 
H 

S'  +  x  sin  X  -  y  c0S  x 
H 

zh  -  i 

X  cos  X  +  y  sin  X 
-  y  cos  x  +  x  sin  X 
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.  V* 


y2  ^  52 


2  +  ZL  +  z 
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27 

28 

29 

30 

31 


M  =  (xX  +yV  +  zZ)/s  aQ 
s 

M  =  (xX  +9Y  +  iZ)/Sa0 
s 

o>  F, 

1  n 

b\  F  =  -  kfl  a  ^sink,-^  ,Fi)c“k* 
b)  i  kn  n 

a)  f  ,  F.COS0  +  F  sin  *'^"8^,  cos*' sin  8 

b)  F  =  F  «  cos  <$>'  -  Fji  s'n  $ 

,  /  =  ."f  ,  sin  e  *  F  sin  *' cos  8  +  Fz,  cos  *' “»  0 

c)  rz  x  y 

„)  F  =.  F^,  cos  0  -  Fx,  8  sin  0  -  (F^*'  *  %'>  sin  ♦  S'"  ® 

X  +  /z,  +  Fy,V)  cos*'  sine  +  Fy,8  sin  V  cos8+  F^Scos  ♦  cos  8 

b)  Fy=  (Fy,-F>')«sV-V'  +  Vs^' 

c)  l  .  -F  ,sin0-Fx,8cos8  +  (Fy,-Fz,V)-t'^e 

2  +  (F  !v  +  F  .)  cos  *' cos  6  -  Fy,6sin*' sin8  -  Fz,8cos*‘ sin  8 
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a) 

b) 

c) 

a) 

b) 

c) 


F  cosx+F  sin  X 
x  y 

F  sin  X  "  F  cOS  X 

x  y 

-  F 

z 

• 

F  cos  X  +  F  sin  x 
x  y 

F  sin  X  -  K  cos  x 

x  y 

• 

-  F  . 

z 
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34 


CT  =  M  +  o((l  -  M2)/* 


X*  *  XT  *  XA 
s 


xm  ■  xt  '  Vy 

S 


Zm  =  XT  -  Z/z 


Ap  = 


■M  T 
s 


*  (FX+  FXXm)+?  (FV+  FY  •  YM  )  +  £  (FZ  +  FZ  ZM 


4ir  (1  -  M)2  a  S2 
s  0 


Pj(0  =  p.(0  +  Ap 


a)  C  =  2B 

Pjl< 


p.(m)  cos  (m  -  1)  k  •  2  it  B/M 


'  > 

b)  S  =  2B  ^  ^  p,  (m)  sin  (m-1)  k  2irB/M 

Pjk  J 


c)  p..  =  /  C2  + 

Jk  V  Pjk 


pjk  =  20  ,09io  pjk  +  124-58- 


fk  =  k  iiB/2ir 
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SEQUENCE, 

peratqr 


16 


program  HERON  1 

• 

COMMON  Kl,K2,KSr,ASX,Dr'SI,M,N 

COMMON  BET(72>iZET(72)iFN(72«12«#FT(72#12)  .DMDI72)  ,D(22> 

COMMON  PS  1 0 1 K  i  C.OM  ,  AH 

COMMON  XOi YO,ZO,XDH,YDH#ZDH 

COMMON  XHO«  YHO,ZHC,THO,THO,PH!OiPHID 

COMMON  AO  (  72  ) , AK ( 36. 72 > ,  BK (  36 , 72 ) , f ( 3  ) , FD ( 3 > 

COMMON  XD ( 72 )  i  YD ( 72 ) , ZO ( 72 > .  HP < 3 > ,  HOP ( i ) 

COMMON  WDDPI3) 

COMMON  BETP(72,12)«ZETP(72,12> 

DIMENSION  CARD ( 3Q ) , RN ( 12 ) 

DIMENSION  FF ( 72  > , CPC ( 72  >  ,  SP J ( 72  > , P J( 72 ) , ALPHAJ( 72 ) 

C.PJM (72) 

EQUIVALENCE  (  T  AUO.  T  AUO  ) »  (  XD,  CXD  ) ,  (  Y  Lr ,  CVD  >  •  I  ZD.  CZD  > ,  <  D  <  1  > ,  RN  ( 1 )  ) 
EQUIVALENCE  (CYDH.YDH),  (CZDM.ZDH) ,  ( BET#  rn»  (ZET#CPJ) 

e 

PI«3 ,1415926936 
WR I T  E ( 6 1 , 1  ) 

RE  AC  (  6  0  i  2  )  C  ARD 
WR  J  TE  <  61 »  3 ) C ARD 
RE  AC (60,1027)  I H, JCON,C,RHO 
DO  COO  II. 1,  1H 

RE AC  <6OilO28)OMN#R,0B,EPF(EPL,PSIO 
0M«CMN*6, 203185307/60, 

REAC (6O1IOOO  )M,N  , I ND 1  KSFtKSX  ,K1,K2 

NX«MAX1(KSX,K1iK2>*1 

kf  »ksf*i 

AM«XM«FLOATF(M) 

ZZ«2 ,/XM 

DT  .6,283105307/IOM‘AM) 

IT  .0,2831053O7/(OM*B0) 

DPSI.6.283185307/XM 

TT.TT-DT 

MFT.M-1 

DMDUJ.O, 

DDSI»360,/XM 
DO  409  J.l.MPT 
405  DMD( J*l>«DMn<J)*DDSI 

REAC  (60,1028)  <D(n,|«l,N) 

00  10  I  * 1 , N 

10  RM  I  )  «R*D(  I  ) 

REAC  <60it028)XH,rH,ZH,TH,PHI,XDHl YDH.ZDH,THD,PHID 
IF<  XDH.EQ, YDH, AND.YOM.EQ.ZDH, AND.ZDH.EQ.O, >1083# 1064 
1 0  83  V  «  0  . 

QO  TO  1085 

|084  V«SCRTIXDH*XDH«  YDH*YDH*ZOH*ZDH > 

1089  CONTINUE 

I  9* I F I  X ( B8 ) 

WRITE (61, 4)  I9iR,0MN,EPF,EPL.PS1O,TH,PMI ,XDH, YDH,ZDH,  V 

WRITE(61,5)XH,YH,ZM 

MR  I TE ( 61 , 6 ) 

SINPD.SINF(PHID) 

COSP  D«C0S ( PH  I D ) 

XHO»XW 
YHO« YM 
ZHO»ZH 
TH0»TM 
PHIC.PHI 


READ  1 
READ  2 
READ 


READ  4 

READ  5 


37 


CHJ«ATAN(YOM/XDH) 

irtKDM.KO.O.Illill 

11  ir(YDH,BQ,0, >14,15 

14  CHI*0, 

00  TC  12 

15  CWJiPI/2, 

12  COSCH! ■COS(CHI ) 

S!NCH!»$IN(CHI ) 

HR | TE ( 61, 9)N 
NR|TE(61.61)M 

8  N  IS  NUMBER  Of  RADIAL  ST  AT  JONS 

8  M  is  the  number  or  azimuth  stations 

8  X,  V,  Z  IS  THE  POSITION  OP  OBSERVER  ********* 

8*  bet  is  the  flapping  angle,  ******* 

8*  ZFT  IS  the  lag  angle,  ******* 

8  BETA  AND  ZETA  ARE  INPUT  IN  DEGHEES, 

GO  TO  (16,18)  IND 
16  CONTINUE 
HR  I TE ( 61 , 7 ) 

READ  (60,1028) ( BE T  (  I  )  ,  ] » 1 , M )  READ  6 

READ  (  60,1028  ) ( ZET  (  I  >,  I *1, M)  READ  7 

HR  I TE ( 61 , 67 ) 

WRITE (61, 68) ( DMD (  I  ),8ET< l ),ZET{  I  ),  !«i,M) 

GO  TO  20 
18  CONTINUE 
HR  I TE ( 61 , 8 ) 

CALL  INPUT  (1) 

20  CONTINUE 

DO  25  I  3*1  *  N 

BET(I3)*BET(I3)*0, 01765329 
25  ZET( I3)»ZET(I3)*0, 01745329 
REMIND  2 

GO  TO  (30,40)  IND 
30  DO  35  J*1,N 

9 

8  The  FOLLOHING  STATEMENTS  INPUT  THE  ETAN  AND  ETAT  ARRAYS  JNTO  THE  FN  AND 

9  FT  STORAGE  LOCATIONS,  THESE  SAME  LOCATIONS  WILL  LATER  BE  USED  FOR  THE  FN 

8  AND  FT  ARRAYS, 

8 

READ  (60,1028)  ( FN ( I , J  )  ,  I  «1 ,  M  )  READ  8 

35  READ  (60,1028)  ( FT ( I , J ) ,  I  *1 , M  )  READ  9 

HR  I TE ( 61, 69 ) 

HR  I TE ( 61 , 63 ) 

HR  I T  E  <  61 , 64  > (RN(  I  ),  !»1,N) 

DO  111  J«1,M 

HR  I TE ( 61 , 65 ) ( ( DMD ( J) , ( FT ( J, I ) , I *i*  N)  > ) 

111  HR  I TE  <  61 , 1091 ) 

HRJTE(61,71) 

HR  I TE ( 61 , 63 ) 

HRITE(61,64) (RN(  I  >,  1*1, N) 

DO  112  J«1,M 

HRITE(61,65)((DMD(  (FN(J,I), J«1,N)  >> 

112  WRITE (61*1091 > 

GO  TO  42 

40  CALL  INPUT  (3) 

42  DO  70  15«1,N 
PS  I  • 0 , 0 
DO  <0  I4*l, m 
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COSB<COS<06T< !4) ) 

C0SPZ<C0S(PSl*ZET(|4) ) 

SINFZ<SIN<PSI<ZET< 14) ) 

S I NG  <S I N (BET ( 1 4 ) ) 

COSP  «C0SF ( PS  I ) 

S I NP  »S l NF  <  PS  I ) 

8 

«•  RELATIVE  COORDINATE  OF  BLADE  STATION,  • 

8 

IF (EPF ,0T,EPL)54,58 
54  T«<RN<  I5).£PF)*C0S8*EPF.EPL 

XD(  I4X»EPL*C0SP-T*C0SPZ*FN<  14,  I5>*SJNB*C0SPZ 
C*FT  ( !  4 ,  1 5  )  *S 1 NPZ 

YD(|4)t»EPL*SJNP«T*SJNPZ»FT<I4, I  5 ) <COSPZ<FN<  1 4 ,  1 5 ) *S l NB<S I NPZ 
ZD<  |4X>(RN<  I5)»EPF><SJNB-FN(  14,  I  5  ) *COSB 
GO  TO  60 

58  XD< I4>..<EPF*GPL<0QSB)*C0SP.<fiN< I  5 > *EPL > <CQSB<COSPZ<FT <  t4, 15  I* 
1S1NPZ*FN<  14, J5)*S!N8*C0SPZ 

YD( 14).*  <EPF<EPL<COS0 )<SJNP.(RN< I  5 > <EPL ) *COSB*S I NPZ -FT < 1 4 , I  5  >  *COSP 
1Z*FN( 14, J5)*SINB*SJNPZ 
ZD(  I4)f(RN(|5).EPF(<SINB-FN(  14,  I5)<C0SB 
60  PSI«PSI<CPSI 

WRITE  (2,1003)  ( XD(  I  4 ) , I 4« 1, H ) 

WRITE  (2,1003)  <VD(I4), I 4«1 , H> 

70  WRITE  (2,1003)  <ZO<  14), 14<1,M> 

DO  93  ! 5<1, N 
DO  93  I4<1,M 

• 

•<  AERODYNAMICS  force  component,  <*♦.***.***  < 


f  XOiFXC  ,  YD»FYD  ,  ZD<FZD  IN  THE  FOLLOWING  COMPUTATIONS, 

8 

IF  (N.EQ.l)  91,83 
83  I F ( I5.E0, 1)84,88 

84  DNDR»(FN< I4,2)-FN<  14,1))/(RN(2)-RN(1)) 
DTDMtFT(!4,2>*FT(I4,m/<RN(?)«RN(l)> 

GO  TO  92 

88  I F  < 1 5 , EQ , N ) 89, 90 

89  DNDR«(FN( 14, 15).FN(  14,  J9-1) )/(RN( I5).RN< I5<‘ ) > 

DTDR<(FT< 14, 15 ) -FT «  14,  15-1) )/<RN( I5).RN( 15-1) ) 

GO  TO  92 

90  DNDR<(FN(  14, 15<1).FN( 14, 15*1)>/<RN( I5<1)<RNI I5<1) ) 

DTDR ■ ( FT  ( 14, I5<1)-FT( 14, 15.1))/ (RN( 15*1) -RN( 15-1) ) 

GO  TO  92 

*1  DNDR  <DTDR*0 , 

92  BETP(  14, 15)<8ET( 1 4 ) < AT AN ( DNDR > 

93  ZETP(  14, 15)*ZET(|4)<ATAN(DTDR> 

GO  TO  (94,96)  I  NO 

94  DO  95  J<1,N 

READ  (60,1028)  <FN< I , J) , J <1,M) 

95  READ  (60,1028)  ( FT ( I , J ) , I <1, M > 

WR I TE ( 61 , 62 ) 

WR l TE (61 » 63 ) 

WRJTE<61,64)(RN<I),  l <1 , N ) 

D01112  J»1,M 
WRITE  (61,1091) 

1112  WRITE(61,65)((CKD(J),(FN(J,I),I<1,N)  )) 

WR|TE(61,66) 

WR I TE ( 61 , 63 ) 

WR|TE(61,64)(RN(|),  I  <1 , N ) 

DOllll  J ■ 1 » m 
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WR 1 T E ( 6 1 , 1 0  91 ) 

till  HRITE<61.65> MOMDf J>» <FT(J, I), 1-1, N>  )) 

GO  TO  97 

96  CALL  INPUT  (2) 

97  00  160  15-1, N 
PS  I  *  0  . 

DO  100  14-1, N 
ZE-7ETPI 14, 1 5 ) 

B£»5ETP( I4| 15) 

XD(!4)»FT(I4,m*SlN(PSI*ZZ  )*PN(  14,  (  5  )  -S I N ( BE  )*COS(PSI*ZE  ) 
YD  ( 1 4 )  «*FT ( Mi J5)*C0S(PSJ*Z6  >*FN(I 4, Z5>*SIN<Bfc  >*SIN»PSJ*ZE  ) 
Z0( I  4  » ■ -FN ( 14, I5)»COS(0E  ) 

100  PSI-PSl*DPSf 

WRITE (2. 1003) <XD( J4>, J4-1.M) 

WRITE (2, 1003)  <YD( 14) , |4«1, M> 

WRITE (2. 10 03)  ( ZO ( I  4 ) , |4-1,M) 

160  CONTINUE 

END  FILE  2 
REWIND  2 
9 

•*  CALCULATE  HARMONICS  •  •••***•*..♦*♦** 

C 

17*0 

DO  190  12-1,2 
GO  TO  (165,166)  12 
1*5  KSPaKX 

GO  TO  167 
166  KSP-KF 

1*7  IF  (KSP.LE.O)  168,1*9 
1*6  KSP»1 
1*9  DO  190  13-1, N 
00  190  14-1,3 
17-17*1 

READ  (2,1003)  (BET ( 15), 15-1, M) 

SUM1-C, 

DO  172  19-1, M 
172  SUM1-SUM1*BET< J9) 

A0(I7)-SUM1/XM 
DO  190  16-1, KSP 
XX-CPSJ*FL0ATF(I6) 

SUM2-SUH3-X-0, 

DO  180  18-1. M 

SUM2-SUH2*BFT( 18)«COS(X*XX) 

SUM3-SUW',*8ET(  I8)-SIN(X*XX) 

8  WRITe(61,1044JSUM2,SUM3 
ISO  X-X«l, 

AK( 16, J7)-SUM2*ZZ 
BK( I  6 , I  7 ) -SUM3-ZZ 
190  CONTINUE 
I 

••*«**.  ********  BEGIN  PART  2  ---  CALC,  OF  SOUND  FIELD  ••***•*•***** 

« 

1 3- I F I X ( BB ) 

HR  I TE ( 61 , 78 ) 

DO  700  1  1 1 » 1 , JCON 
17.0 

REAP  (60,101*)XOP, YOP.ZOP 
SXBW«*XOP*XHO 
SYBW-- Y0P*YH  0 
SZBW--ZOP*ZH0 

IF (IXBH.EO.SYBH, AND, S YBH , 10 , SZBW , AND, SZBH.EO, 0, ) 1073,1072 
1071  SW-0. 
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00  TO  1075 

1079  SM*SQnT<SXBH*SXBH*SYBH*SYBH*SZBH*SZBH) 
1075  00a  SXBM»XDM«SYBH*YDH*SZ0H*ZDH 

QTaC*C«V*V 
QS»CQ*QQ*SH»SH*QT 
I T  <  CS , EO i 0 i >1071,1092 
1071  TAUCPa«QO/QT 
GO  TO  1069 

1092  T*L0P»<-00*S0RT(0S> >/OT 

1069  CONTINUE 

XPH*XNO-XDH*TAUOP 
YPH« YNO-YDH*T*UOP 

zpm»zno-zdh«tauop 

3X0PH»$X8H*XDH*TAUOP 

SY0PH*SY8H*YDH*TAUQP 

SZ9PH*SZBH*ZDH«T  AUOP 

BHO«SXBPH*xnH*SYBPH*YDH*SZBPH*ZnM 

BM0»9P0/C/C/TAU0P 

TTa« ,2831853/(0M.8B> 

DT*6,2831853/(0M*4M) 

T*0 , 

MBaX/JF IX(BB) 

DO  A  0  0  MSel.MB 
I  7a  I  7 ♦ 1 

XO«XOP*XDH* ( T  -  T  T / 2 , ) 

YO»YOP*VDH. ( T-TT/2, > 

ZO*ZOP*ZDH* ( T-TT/2, ) 

PJ(  I7)«0, 

DO  350  I2«l.  13 

PN6uaPSI0*6,28316t  •  F L 0  A TH I  2  *  1 > / R B 
T  AU1 ■ T  AUOP 
11*0 
JJ«3*N 

DO  ’3t  M»i,N 
C 

6*  ,TE  RETARDED  TIME  *  •  *  * 

8 

-  J 

i\K  *  1 

210  TAUC*TAU1 
TPaT.TAUO 

TAUlaSFITP.KK,  I4,PNEB,CHl )/C 
1*1*1 

IF  (I.GT.50)  260,220 
220  ERaABSF  < T A  i.TAUO) 

IF  ? ER , GT  ,  ,  0 01  )  225,228 
225  KK»1 

GO  TO  250 

228  IF  (EB.GT,  ,  0001  )  230,232 
230  KK a 3 

GO  TO  250 

232  IF  (ER.GT,  ,  00005  )  234,236 
234  K K * 5,  0 

GO  TO  250 

236  IF  (ER.GT,  ,00001  )  238,242 
238  K K  s 2  0 

GO  TO  250 

242  IF  (KK.EQ.KX)  260,254 

250  IF<kx«KK)251,210,210 

251  KKaKXf  GO  TO  210 
254  KKiHXl  GO  TO  210 

360  CONTINUE 
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T  P ■ T »TAU1 

8 

8  BEG I N  CALCULATION  OF  VELOCITIES  AND  ACCELERATIONS, 

C 

PSI»PN£W*0M*TP 

XN«XH0*XUK«TP 

YM*YHC*YDH»TP 

ZH«ZH0*ZDH*1? 

TH*TH0*THD*TP 
PHl«PNIO*PHID*TP 
COST  *C0SF  <  T W  > 

S I NT*S  J  NF ( TH ) 

TANP*TANF<PHI) 

PHIF*ATAN(C0ST*TANP> 

SINPHIP,S1N<PHJP) 

COSPHlP»COS<PHIP) 

PHICP*COSPH!P*COSPHIP*»PH!D*COST/(COSF(PHI )**2)-ThD*SINT«TANP> 

10 AS  F ORMT ( 6X ,  6HPS I X8 « , 12E10 , 5 > 

DO  265  15*1,3 
11*11*1 
WP(  !5>«AO<  ID 
X  K  » 1  , 

WDP ( 15) * wDDP  < I  5 ) *  0 , 

DO  263  KL'l.KX 
ARO*XK*PSI 
COSHP«COSI ARQ) 

S I NKP»S I N ( ARQ ) 

HP  ( I  5 ) * WP ( l  5 ) * AK ( KL i  I  I ) *C0SKP*bK ( KL ,  I  I  )*‘  >P 

HOP! 15)* WDP ( I5)-XK*0M*( AKIKL. ! I >  *S  J  NKP  •  RK  UL , I  I  )  *COSKP) 

WDDP( 15) *WDDP( 15 ) -XK*XK*OM*OM*(AK(KL, I  I >*COSKP*BK(KL,  I  I  )*SINKP) 

263  XK.XK*1, 

265  CONTINUE 

• 

THD2»THD*THO 
P2*PH I DP*PH I  DP 

X*WP(1 )*C0ST*WP(2)*SINPHIP«SINT*WPI3)*C0SPHIP*SINT 
Y*WP(2)*COSPHIP-WP<3)*SINPHIP 

Z**WP(1 ) *SINT*WP<2) *SINPHIP*C0ST*WP( 3) *COSPHlP*COST 

0 

XD*NDP(1)*C0ST»WP(1)*THD*SINT*(WP(3)*PHIDP«WDP<2)  )  * S  I  NPW I P*S I  NT* 
1(WDP(3)*WP(2)*PHIUP)*C0SPHIP*SINT*WP(?)*THD*SINPHIP«C0ST*WP(3)* 
2THD‘C0SPHIP*C0ST 

YD«(WDP(2)-WP(3)*PHIDP)*C0SPH!P*’(WP(2)*PHIDP*WDP(3)  )  *  S  I NP  W  I  P 
ZD*.WCP(1)*SINT-HP<1)*THO*COST*{WDP(2>-WP(3)*PHIDP)*SINPHIP*COST* 
1(WP(2)*PHIDP*WDP(3)  )*C0SPHIP*C0ST-WP(?)*THD.SINPHIP*SINT.WP(3)* 
2C0SPH I P*S I  NT 

0 

XDD*(WDCP(l)=WP(l)*THD2  )*C0ST»2,*rt0P(l)*THD*SlNT*(2,*WP(?)* 
1PWICP*TND*2,*WDP(3)*THU)»C0SPH1P*C0ST*(WDDP(2)®WP(2)*(P2* 

2TMQ?  ).2,*WDP(3)*PHIDP)«SINPH!P*SINT*(21*WDP(2)*PNI0P*WDUP(3)- 
3WP (3 ) * (P2  ♦ THD2  ) )*C0SPHJP.SINT*(2,*W0P(2)»THD-?,*WP  (3) 

a*phidp*tho)*sinphip*cost 

YD0«(WDCP(2)*wP(2)*P2  )*C0SPHIP-(WUDP(3)-WP(3)*P2)*SINPMIP 

ZDD»-  2,*WUP(1)*TH0*C0ST=<WDDP<1>-WP(1)*THD2  ) *S I  NT* ( 2 , *rfDP ( 2 ) 

l»PHIDP*RDDP(3)eWP(3)*(P2*THD2  ) ) *COSPH I P* S I  NT .  (  2  ,  * WDP < 2 

2)*Tt-D®2,*WP(3)*PHIUP*THD)«SINPHIP*SlNT»(wP(2)*PHinP*WnP(3))*COSPHI 
3P«SINT*2,*THD*(WDDP<2)*WP(2)«(P2*THD2  > -? , *WDP  (  3 ) *PH  I  DP 

A  )  *S I NPH  I P»COST 

a 

XRARS  XO-XH-X*COSCHl *Y*S INCHI 
YBAP*YO-YH-X*SINCMl*Y*COSCHl 
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ZBAR*ZO*ZH*Z 

a 

CXO»XCH*XD*COSCHJ*YD*SlNCHI 
C YD*  YDM*XD*S!NCH!*VD*COSCH! 

C7D*  ZDN*ZD 

a 

CXDC«XDD*C0SCH1  *YDU*SINCH! 

CYDC»-YDD*COSCH!«XDD*SINCHI 

CZDC«*ZDD 

a 

SfS*SQfiT(XfiAR»XUAt..rHAR*YBAR*ZBAB*ZaAfl) 

S 

XNA5*»CXD*CXD*CYD*CyD*CZ0«CZD)/C/C 
XHS*(X«AR*CXD*XBa»*CYD*ZSA«*CZOI/(SS*C) 
XMD««<X8AR*CXD0*YBAR*CYDD*ZBAR*CZDDJ/f SS*C) 

0 

0  BEGIN  CALCULATION  OF  FORCES  ANU  DERIVATIVES, 

0 

DO  275  15*1,3 

JJ»J J*1 

F(  1 5  )  ■ AO ( JJ) 

F  D  (  I  5  >  ■  0  . 

XK«  1  , 

DO  270  16*1, KF 
COSK*COSF(XK*PSI ) 

SINK*SINF(XK*PSI ) 

F< 15)»F(I5)*AK(  16, JJ)*COSK*BK< 16, JJ>*S|NK 
FD(15)«FD<I5).0M«XK*(AK(I6,JJ)*SINK-BK(I6,JJ)»C0SK) 

270  X  X ■ X  K  *  1 , 

279  CONTINUE 

0 

FX*F(1)«COST*F<2>*S!NPH!P»S!NT*F(3)*COSPHIP*SINT 

FY*F(2)*C0SPHIP-F<3)*$INPWJP 

FZ»-F(l)*SlNT*F(2)*SINPHIP*COST*Ft3)*COSPMlP*COST 

0 

FDX*FD(l)*COST.F (1 J*TWU*SINT.<F  < 3 > *PH I  DP -FD < 2 > )*SINPHIP  *SInT* 

1  <FD<3)*F<2)*PHIDP)*C0SPH!P*S JNT*F(2)*THD*SINPHIP* 

2  COST*F(3)*COSPH|P«COST«THD 
FnY«(FD<2)-F(3)*PHl0P)*C0SPHIP-(F(2l*PHIDP*FD(3) >*SINPH!P 
FOZ**FD<l>*S!NT"FU>*THD*COST*<FD(2>«Ft3>*PHIDP>*SlNPHlP*COST* 

1<F<2)«PMIDP*FD(3) )*C0SPHIP«C0ST*F{2)«THD*SINPMIP*SINT.F(3)*C0SPMIP 
2*SINT*THD 

8 

FX*FX«COSCHI*FY*SINCHI 

Fy*FX*SINCH]»FY*COSCHJ 

FZ*.FZ 

0 

fdx«fcx*coschi *FDY*S INCHI 

FDY«FCX*SINCHl.FDY*COSCNI 

FOZ.-FDZ 

0 

0  BEGIN  SOUND  PRESSURE  CALCULATIONS, 

XT*XMDS*C«(l,‘yMAS)/SS 

XHXS*XT-XD/X8AR 

XMYS*XT-YD/Y8AR 

XmZS*xT«ZD/ZBAR 

9 

TFMP  *1 ,  -XMS 

TEHP2*TEMP*T6HP 

TMP*-1,/(4,*P|*TEHP2*C*SS*SS) 

DP»(X9AR*UDX*FX*XMXS/TEMP)*YBAR*lFOY*FY*XMYS/TiLMP)*ZBAR*(FDZ*FZ* 

1XMZS/Tr-"PI)»TMP 
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rlNfO« 


PJ( I  7  >  «P J t 1 7 )  *DP 
330  CONTINUE 
350  CONTINUE 
T«T*DT 

400  CONTINUE 

6 

8  HARK  ON  I C  ANALYSIS  OF  SOUND  FIELD, 

8 

Fia2,«BB/XM 
F2a2,*Pl*B8/XM 
F3«CM«89/ ( P  I  ♦ P  I  ) 

F4a20,*, 43429448 

XM0.XH/R8 

BMaXMB/2, 

BM2»FL0AT(MB/2) 

IF  (BH,GT.BH2I  605,610 
605  KMAXbM/<2*IFJX(88)  ) 

QO  TO  615 

610  KMAXaN/(2*IFIX(»B>  >«1 
615  XL«C, 

DO  <80  L1"1,KMAX 
XL**L*1, 

SUMlaSuH2«X*0 , 

DO  <60  HMal.MB 

SUMlaSUHl*PJ<Hh)*COS<XaF2*XL) 

SUM2aSUH2*PJ<HM)*S|N<  X«XL*F2) 

664  FCRMAT(2X,2I2,9(JX,F10,5)> 

660  X«X*1 , 

CPJtLl ) *F1*SUM1 
SPJ(L1 )«F1*SUM2 

ALPNAJ(LHaATAN«SPJ(Ll)/CPJ(Ll)  > 

I F ( CP J ( LI ) , EO , 0 , >661,662 

661  AlPKAJ(Ll)aPI/2, 

662  SPW(L1)«ABSF(SPJ<L1)  ) 

CPHLl)«A8SF(CPj(Ltn 

PjK(Ll)'SORT(CPJ(H)»CPJai)*SPJ(Ll)*SPJ(U)  > 

PJK<L1 )aF4*ALOQ(PJK(Ll» >*124,58 
FF(Ll)aXL*F3  /<l,aBNO) 

660  CONTINUE 

WRITE  (61,1030) 

WRITE  (61,1032)  III, XOP, YOP, ZOP 
WRITE(61,1029)XPH,YPH,2PH 
WRITE  (61,1030) 

WRITE  (61,1033) 

WR I TE ( 61 , 79  ) 

DO  700  Ll'l.KHAX 

700  WRITE  (61,1034)  LI , FF ( LI ) , CPJ< LI ) , SP J{ LI > , PJK ( LI ' 

WRITE  (61,1030) 

800  CONTINUE 

WRITE  (59,1002) 

WRITE  (61,1030) 

WRJTE  (61,1002) 

FORKAT(1H1,60X,16HPROGRAH  .HERON  1) 

FORKATt  1 0  AQ ) 

FORKAT (3(28X,10A8/) ) 

FORKAT(1HO,28X,63HBLADE  LOADING  AND  NOTION  DATA  (PARALLEL  AND  NORM 
1AL  TO  THE  SHAFT//) 

4  FORK  AT  ( 1H  ////,  28X,  22H  GEOMETRY  OF  ROTOR  H  ,,U0  ,11X 

1  ,6W8LA0ES/,IH  ,50X,E10,3,3X,3HFT,,4X,6HRADIUS/,1H  ,  50X, El  0 , 3 , 3X , 
26HR,P,M,/,1H  ,50X,E10,3,3X,3HFT,,4X,17HFLAP  HINGE  QFFSET/.1W  ,50X, 
3E10,0,3X,3HFT,,4X,17HDRAG  HINGE  OFFSET/.1H  , 50X , El  0 , 3, 2X , 7H0E0REES 
4.14W  AZIMUTH  PHASE/, 1H  , 50X, E10 , 3, 9H  DEGREES, 27H  REARWARD  SHAFT  I 


44 


D  o  n  Q  o 


5NCL  I  NAT !0N/,1H  ,50X610.3,98  DEGREES , 318  SHAFT  INCLINATION  TO  STAR 
6BOARD/,50X,ElO,3,2X,6HFT/SfcC,2!jH  VELOCITY  IN  X-D I RECT 1  ON/ 1 1H  ,50X 
6,E10,3,2X,6HFT/SEC,25H  VELOCT I Y  In  Y. DIRECT  ION/, 1H  , 50X, E10 , 3, 2X , 

86HFT /SEC , 25H  VELOCITY  IN  Z-0 I HfcCT 1  ON/ , 1H  , 50 X , ElO , 3 # 2X , 6HFT/SEC,  2 
90H  RESULTANT  VELOCITY/) 

5  FORFATdH  ,50X,E10,3,3X,3HFT,,18H  X  HUB  ORDINATE/, 1H  ,5QX,E10,3 
1 , 3  X , 3NF  T , , IflH  Y  HUB  ORDINATE/, 1H  ,  5 0 X ,  F 1 0 , 3 , 3 X , 3HF T  , , 18H  7  H 

2UB  CRCINATE) 

9  F  ORF  AT  < 1H  »  ?8X , 27HNUM8ER  OF  LOADING  STATIONS*, J6//) 

61  FORFATdH  , 28X, 37HNUMBER  OF  AZIMUTH  INTEGRATION  PO I  NTS* ,  1 6// > 

7  FORFATdH  , 26X , 49HFORM  OF  INPUT,  SPANW I S6/ AZ I MUThW I S6  DISTHUBUTION 
IS//) 

8  FORFATdH  ,  28X ,  23HF0RM  OF  INPUT,  HARMONIC//) 

67  FORFATdH  ,//28X,l3HAZJMUTH  ANGLE , 7X ,  1 0 HF L AP  ANGLE, 10X,  9HLAG  ANGL 
IE//  ) 

68  FORFATdH  ,28X,El0,2,10X,F7,3,13x,f  7,3) 

69  FORFATdH  , ///28X,57HBLADF  ELEMENT  ELASTIC  DISPLACEMENTS  PARALLFL 

i to  the  shaft//) 

63  FORFATdH  ,  2X,  13HAZIMUTH  ANGLE  ,  30X,  1  AHR/.D  I  AL  STATION//) 

64  FORFATdH  ,15X,10F10,3) 

65  FORFATdH  .  5X,  F10.3  ,10610, 3) 

71  FORFATdH  , ///28X.57HBLADE  ELEMENT  ELASTIC  DISPLACEMENTS  NORhAL 
ITo  THE  Shaft//) 

1091  FORFAT(Z) 

62  FORFAT (1H1  j  28X , 26HBLAUE  THRUST  LOADING  (LBS)//) 

66  FORFATdH  , ///28X,  lflHBLAUb  DRAG  LOADING//) 

78  FOOLJddH  .///A8X.20HCOMPUTFO  SOUND  FIELD//) 

79  FOHFATdH  ,  28X  ,  8HHARHON  I  C  ,  6X  ,  9HFRE  QUENC  Y  ,  8  X  ,  1  4  H  I  N  PHASE  HRES,,5X 
1.10FOLT  OF  PHASE  PRFS,  ,  3X, 20 H SOUND  PRFSSURE  LE VE L/ # 3 Q X , 3HN0 , , 31 X  , 

23H(E  ,  S  ,  F , ),10X,8h(P,S,F, ) ,9X,28H(DB,  RE,  ,  0002  OYNES/SU  ,  C« , ) ) 

1029  FORFATdH  , ///28X,  25HRETARUED  HUB  CO. OH D I N A TES , 1 0 X , 4 H XPH* , El  6 , 8 , 1 X  ,  4 H YH 
1, 4H Y PH s, £  16, 8, 1X.4HZPHI,  6 16, 3) 

1032  FORFATdH  , ///28X.21H0BSERVER  CC*OHD I  NATES  »  4X, 

U  2Hj=,l4,5x(Jhxa=,E1618,?x,3HYOs,E16,B,2X,3HZ0=,F16,P) 

1000  FORFAT(20I4) 

1002  FORFAT  (2X.10HEND  OF  RUN) 

1003  FORFAT!  BE  1 6 , 8 ) 

1016  F  ORA  AT  (  BE  1 0 , 0 ) 

1 0?  7  FORFAT(2I4,?X,6F10,0) 

1028  FORAAT(RF10 ,0 ) 

1030  FORMAT  (/) 

1034  FORFAT(30X, H.4E20.8) 

10  33  FORF AT (32X, 1 HK, 13X, 4HF(K ) ,  1 3 X , 6HCPJ C K  >  , J 5 X , 6HSP J ( K  > , 1 4 X  ,  5HP J ( K) ) 

END 

SUBROUTINE  INPUT  (II) 

C  0  H  f  o  K  K1,R?,KSF,KSX,DPSI,M,N 

COMFOK  °ET(/2),ZET(72)IFN(72,12),FT(7?,12)  ,DMD<72)  ,D(2P) 

DIMENSION  ZH0?<12)»CFN(12,lP),/P03d2>,cmi2,l2>,E<l?>,SFNU?.12> 
DIMENSION  ZR01(1),CHET(72)  ,GGGU),SRET  (  7  2  )  ,HHH(1  )  ,  SFT  ( 12, 12) 

call  parameter  or  i  returns  beta  and  zeta  arrays  only, 

Call  PAPAMF'ER  OF  t  RETURNS  FN  AND  FT  ARRAYS  ONLY, 

CALL  PARAMticR  OF  J  RETURNS  ETAN  AnD  R  T  4  T  ARRAYS  ONLY, 

WRITE  (61.1505) 

1  =  0 

77  FORMATdh  ,  ?8  X  ,  3  4HBL  ADE  LAGGING  HARMONICS  (  DFGREES  )  ) 

GO  TO  (50,110,120)  II 
50  CONTINUE 
C 

C  COMFUTE  BETA  AND  ZETA  ARRAYS, 
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non  c  no 


a 

RE  AC  (60,1028)  ( CflkT (  KK ) , KK ■ I »  Kl ) 

XJ*1  , 

74  F0RMAK1H  ,  3  2  X  1 1 H  A  »  I  2 , 7  X  ,  E10,J) 

76  FORMATdH  ,28X,34HBUADb  FLAPPING  HARMONICS  (DEGREES )  > 
WRtTE<61,76> 

WRITE«61,74)(J,CBET( J), J»I,K1I 
RE  AC  (60,1028)  ( SBE T ( KK ) , KK* 1 , Kl ) 

84  FORMAK1H  ,  32X , lHB , 1 2 , 7X ,  E10.3) 

WRITE(61,84)(J,S0ET(J),J*1,K1) 

XJ«2  , 

X«-l  , 

DO  60  11*1, M 
X  »  X  *  1 , 

PS  I  *  X  *  DPS  I 
BEK  11  >«C0fcT(  I  ) 

Y*0  , 

DO  60  12*1, Kl 

Y  *  Y* 1 , 

Y  K  ■  F  S  I  *  V 

60  BFT ( 1 1 ) SBET<  U ) *CRET ( 12) *COS( YK ) *SBFTt 12) *S  !N< YK ) 

completes  computation  or  BETA  ARRAY, 

75  format < ih  , 28x , 3ihbl ade  drag  load  harmonics  ( lbs > / / / j 

05  F  ORM  AT  ( 1H  ,  28X,  33HBLADE  THRUST  LOAD  HARMONICS  (LBS)//'/) 
REAC  (60,1028)  <CRET(J),J»I,K2) 

X  J*3  , 

HR  I TE  <  61 , 7  7  > 

MRITE(61,74)(J,CBET(J),J«I,K2) 

REAC  (60,1028)  (SRET(J) , J»1,K2) 

X  J  *  A  , 

WRITE(61,84)(J,SBFT(J),J*1,K2) 

X«-l  , 

DO  70  I  3  s 1 , M 
X  •  X  ♦  1 , 

PS  I =  X  *  DP S I 
2  E  T ( I3)=CRET( I ) 

Y  *  0  , 

DO  70  KK  s 1 , K2 
Y* Y«1  , 

Y  K  s  Y  *  F  S  I 

70  2FT(I3)*ZET(I3)*CBET(KK)*COS(YK)*SBET(KK)»SIN(YK) 

79  FORMAT  { 1 H  ,32X,1HA,I2,7X,10E10,3) 

78  FORMAT (Ih  ,32X,1hR,I2,7X,10F10,J) 

72  FORMAT (Ih  ,  2 8 X , 1 A  HR  A U  I  A L  S T A T I  ON , F 1 0 , 3 ) 

COMPLETES  COMPUTATION  Of  Zfc  T  A  ARRAY, 

GO  TO  1200 

110  KSTFsKSF 
K  F  s  1 

GO  TO  125 
120  KSTFsKS* 

KF 

125  DO  130  NN*1,N 

REAC  (60,1028)  (CFN(NN, J ) , J* I , KSTP ) 

X  J  g  K 

RFAT  *  (60,1028)  (SFN(NN,J),J=1,KSTP) 

Xo«A  , 

RFAC  (  60,1028)  ( C  F  T (NN,J),Js I ,KSTP) 
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«a  C3  C3 


XJ.7, 

REAE  (  6  0  #  1  0  28  )  (  SF  T  (  NN  <  J )  i  J*  1 ,  i',  S  T  P  ) 

XJ«?  , 

130  CONTINUE 

69  FORMATdH  , ///28X, 57HBLADE  FLEMFNT  ELASTIC  DISPLACEMENTS  PARALLEL 
1  To  the  shaft//) 

71  FORMAT ( 1H  , ///28X.57H0LAOE  ELEMENT  ELASTIC  DISPLACEMENTS  NORMAL 
1TO  THE  SHAFT//) 

I F (  KFjEO.X  ) 1 0 1 / 50  2 

101  WR I TE ( 61 , 85  ) 

GO  TO  503 

502  WR | T  E ( 61 «  71 ) 

503  WBITE<61,72>  (  D(K),Ksl,N) 

DO  104  JsI.KSTP 

104  WfiITE(61,79)(J,CFN(NN,J),NN8l,N) 

DO  105  Jsi.KSTP 

105  WRITE(61,78)(J,SFN(NN,J),NNsl,N) 

I F (  KF , EQ , 1  >102,504 

102  WR I TE ( 61 , 75  ) 

GO  TO  505 

504  WR I T E ( 61 , 69  ) 

505  WR I TE ( 61 , 72 ) (  D(K)«K*1*N) 

DO  107  J  *  I  *  K  S  T  P 

107  WR]TE(61»79)(J,CFT(NN,J),NN81*N) 

DO  108  J  8 1 «  K  S  T  P 

108  WRITE(61,78)(J,SFT(NN,J),NNrl,N> 

126  CCNT1NUF 

DO  140  N  N  8 1 ,  N 
X*0  , 

DO  140  mm«1,M 
PS  I » DPS  I *X 
xax*l , 

FN  (  MM,  NM  *  CFN  (  NN,  I  ) 

FT(mm,NM«CFT(NN,  !  ) 

T  =  0  , 

DU  140  14*1, KSTP 
Y  8  Y  ♦  1 , 

ARG»Y*PS  I 
SlN*RGsSIN( ARG) 

COSARGsCOSI ARG) 

FN(MM,NN)*FN(MM,NN)*CFN(NN,I4)«COSARG*SFN(NN,I4)«31NARG 
140  F  T< MM,  NN ) *FT (MM, NN ) *CFT ( NN  ,  I  4 ) *  COS  A RG ♦ S F T ( NN ,  I  4  )  «S I N A RG 

COMFLETES  FN  and  FT  OR  ETAN  and  ETaT  CALCULATIONS, 

1028  FORMAT  ( 8F 1 0 , 0  ) 

1200  RFTLRN 

1500  FORMAT  ( 10X, 7E16 , 8  ) 

1505  FORMAT  (/) 

END 

function  sf (tp,kk,num,pnew,chi ) 

o 

COMMON  Kl,n2,KSF,KSV,DPSI,M,N 

COMMON  BET ( 72 ) , ZET( 72  )  ,FN( 72, 12) ,  F T  (  72, 1 2 )  ,DmD(72)  ,D(22) 

COMMON  PSI ) , K ,  C,OM  ,AM 

common  xo, ro, zo, xdh, yoh, /oh 

COMMON  XHO,YHO,ZHO,rHO,THD,PHIO,PHlD 
COMMON  A O  (  '2  > , AK (  36, 72) ,RK<  36, 72  )  ,  F  (  3 ) , FD( 3 ) 

DIMENSION  Wi 3) 

0*  CALCULATED  sound  source  to  observer  distance, 

0 
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PS  I *PKfcta*OH*TP 
XH«XHO*XUH«TP 
Y W*  Y  HO  *  Y  DH«  TP 
2H«2HO*?DH«TP 

th»tho*thd«tp 
PHJiPP10*PH!0«TP 
PW1F«ATAN(C0SF(TH)»TANF(PH| ) ) 

DO  100  !  1  *  1  >  3 
I  !»3*(NOM-l)*U 
WP*  A  0  ( I ! ) 

XK»1  , 

DO  80  12*1, KK 

900  F  OR*  AT  (9X.8HSUB  DAT  A , 5X , 7fc 1 4 , 6 > 

WP«uP*AK(|2,m«C0SF(XK*PSI>  ♦BM  I  2  »  !  I  )  *S  I N  (  XK  «PS  I  ) 
BO  XK  *  X  K* 1  i 

loo  w(M)«mp 

COST  «COSF  (  Tm ) 

S!NFHIP»S!NF(PH!P5 

s  tn*»s  I  NF  <  T  H ) 

COSFHlPsCOS(PHIP) 

S ! NP  «S I NF ( CH I ) 

COSP  «C  OSF (CHI ) 

X«H(1)*C0ST*H(2)*SINPHIP*S1NT*w(3)*C0SPH|P*SIxjT 
Y«W(2)*COSPHIP*W(.Y>  •SJF'PHIP 

Z  *  *v>  (1  )«SINT*w(2)»SINPHIP»COST*w(3)«COSPH!P*COST 
Ml ) *XO-XH.X»COSP*Y«S?NP 
W<2 J «YO-YH-X«SINP*Y*COSP 
M3)»?0-ZH*7 

SF»SQPT  <w(1)*W(1)*h<2)*W(2)*h(3)»M3)  ) 

RFURK 

tND 
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TABLE  I. 

MAJOR  FORTRAN  SYMBOLS  USED  IN  PROGRAM 

HERON  1 

(Other  Symbols  Whic 

h  Appear  are  Generally  Used  for  Local  Calculations  and  have 
Self-Evident  Meanings) 

PROGRAM 

SYMBOLS 

ALGEBRAIC  neciMmnM 

SYMBOLS  DEFINITION 

IH 

H 

Number  of  rotors 

JCON 

J 

Number  of  observer  positions 

C 

C'°0 

Atmospheric  speed  of  sound  (ft/sec) 

RHO 

Po 

Atmospheric  density  (slugs/ft3) 

OHM 

°N 

Rotor  speed  (rpm) 

R 

R 

Rotor  radius  (ft) 

BB 

B 

Number  of  blades 

EPF 

6f 

Flapping  hinge  offset  (ft) 

EPL 

*1 

Lagging  hinge  offset  (ft) 

PSIO 

* o 

Azimuth  reference  angle  (radians) 

OM 

0 

Angular  velocity  of  rotor  (radians/sec) 

M 

M 

Number  of  azimuth  stations 

N 

N 

Number  of  radial  stations 

1ND 

IND 

Indicator:  If  IND  =  1,  loading/motion  read  as  span- 

wise/azimuthwise  distributions 

— 

If  IND  =  2,  loading/motion  read  as  Fourier 
coefficients 

hSF 

Number  of  force  harmonics  read 

KSX 

K* 

X 

Number  of  displacement  harmonics  read 
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TABLE  I  -  Continued 


PROGRAM 

SYMBOLS 

ALGEBRAIC 

SYMBOLS 

DEFINITION 

K  1 

Number  of  flapping  angle  harmonics  read 

K  2 

k2 

Number  of  lagging  angle  harmonics  read 

KX 

K 

X 

Number  of  displacement  harmonics  calculated 

KF 

Kf 

Number  of  force  harmonics  calculated 

DT 

At 

Time  increment  between  sound  pressure  values  in  final 
time  history  (sec) 

TT 

T 

1 

Fundamental  blade  passage  period  (sec) 

DP  SI 

Aip 

Azimuth  increment  (radians) 

D(l) 

x(n) 

if^  radial  station  (nondimensional) 

R(I) 

r(n) 

Radial  distance  of  i^  station  -  R  •  D(I) 

XHO 

\ 

) 

YHO 

Yu 

1  Hub  coordinates  relative  to  fixed  ground  axes  (ft) 

Ho 

/  at  time  t  =  0 

ZHO 

ZHo 

I 

THO 

90 

Pitch  angle  relative  to  aircraft  trajectory  coordinates 
(radians)  at  time  t  =  0 

PHIO 

^0 

Roll  angle  relative  to  aircraft  trajectory  coordinates 
(radians)  at  time  t  =  L 

XDH 

*h 

) 

YDH 

Y, , 

1  Hub  velocity  components  in  fixed  coordinate 

H 

/  directions  (ft/sec) 

ZDH 

) 
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TABLE  I  -  Continued 


PROGRAM  ALGEBRAIC 
SYMBOLS  SYMBOLS 


THD 

0 

PHID 

4> 

V 

V 

CHI 

X 

XH 

XH 

YH 

yh 

ZH 

ZH 

BET 

P 

ZET 

5 

FN 

FN°r  nl 

FT 

Fj  or  nT 

XD 

x' 

YD 

y' 

ZD 

z1 

DNDR 

dnN/dr 

DTDR 

dr)j/dr 

AO 

a 

0 

AK 

°k 

BK 

bk 

DEFINITION 


Pitch  and  roll  rates  (radians/ sec) 


Resultant  hub  velocity  (ft/sec) 

Azimuth  angle  between  x  and  X  axes  (radians) 


Rotor  hub  coordinates  at  retarded  time  t'  (ft) 


Flap  angle  (radians) 

Lag  angle  (radians) 

Normal  section  loading  (lb)  or  elastic  displacement  (ft) 
In-plane  section  loading  (lb)  or  elastic  displacement  (ft) 


Coordinates  in  shaft  axis  system  (ft) .  [Also  used  for 
blade  section  force  components  in  same  system  (lb)  ] 


Local  blade  slope  with  respect  to  plane  of  rotation 
(radians) 

Local  blade  slope  with  respect  to  nominal  blade  azimuth 
(radians) 


Fourier  coefficients  of  blade  section  force  and  dis¬ 
placement  with  respect  to  rotor  axes  (ft) 
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TABLE  I  -  Continued 


PROGRAM 

SYMBOLS 

ALGEBRAIC 

SYMBOLS 

DEFINITION 

XOP 

*;<i) 

> 

1 

YOP 

y;(j) 

1 

[  Nominal  observer  position  in  fixed  coordinate 
f  system  (ft) 

ZOP 

ZJ(i) 

> 

1 

SXBH 

*H 

N 

| 

SYBH 

| 

f  Position  of  rotor  hub  relative  to  observer  at  time  t 
>  in  fixed  coordinate  system  (ft) 

S2BH 

N  1 
X 

J 

SH 

SH 

Distance  between  hub  and  observer  at  time  t  (ft) 

TAUOP 

Time  of  propagation  of  sound  from  the  rotor  hub  which 
reaches  observer  at  time  t  (sec) 

XPH 

XH 

> 

YPH 

yh 

1 

1  Hub  coordinates  in  fixed  axis  system  at  retarded 
r  time  t'  (ft) 

2PH 

ZH 

> 

1 

SXBPH 

XH 

SYBFri 

1 

l 

Position  of  hub  relative  to  observer  at  retarded 
f  time  t'  (ft) 

SZBPH 

fH 

J 

1 

BMO 

M 

0 

Hub  Mach  number  component  in  direction  of  observer 

T 

t 

"Observer"  time  t 

XO 

xo 

) 

1 

YO 

Yo 

( 

i 

>  Coordinates  of  "moving  observer"  at  time  t 

ZO 

z 

0 

J 

1 
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TABLE  I  -  Continuod 


TABLE  I  -  Continued 

PROGRAM 

ALGEB'IAIC 

DEFINITION 

SYMBOLS 

SYMBOLS 

2D 

z 

XDD 

M 

X 

YDD 

y 

ZDD 

z 

XBAR 

X 

YBAR 

y 

ZBAR 

z 

CXD 

X 

CYD 

• 

Y 

CZD 

Z 

CXDD 

•  • 

X 

CYDD 

•• 

Y 

CZDD 

•  • 

Z 

SS 

S 

XMAS 

M 

XMS 

M 

s 

XMDS 

M 

s 

F(I) 

F. 

i 

FD(D 

F. 

Aircraft  flight  path  velocity  component  (ft/sec) 


Aircraft  flight  path  acceleration  components  (ft/sec2) 


Coordinates  of  blade  element  relative  to  observer  at 
retc.Jed  time  (ft) 


Velocity  components  of  blade  element  at  warded 
time  (fixed  axes)(ft/sec) 


^Acceleration  components  of  blade  element  at 
retarded  time  (ft/sec2) 


Distance  traveled  by  sound  (ft) 

Square  of  absolute  Mach  number  of  blade  element 

Mach  number  of  element  toward  observer 

Rate  of  change  of  Mach  number  toward  observer 

Aerodynamic  force  component  relative  to  rotor  axes  (lb) 

Rate  of  change  of  aerodynamic  force  component  relative 
to  rotor  axes  (Ib/sec) 
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5.0  PROGRAM  INPUT/OUTPUT 


Table  II  describes  the  preparation  of  the  input  data  cards  for  HERON  1 .  This  is 
fol'^ved  by  an  example  comprising  a  complete  set  of  data  written  on  coding  forms. 
In  conclusion,  the  computer  output  obtained  for  the  example  is  presented. 


TABLE  II.  PROGRAM  HERON  1  .  DATA  INPUT 


T-  r  “ 


u 

a 

o 

o 

O 

CSC 

00 

00 

CM 

CD 

< 

U 

H 

4  4 

1 
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CM 
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CM 
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'■mhU.ILU.U.U.U.U.U.mhmm 


.  .  i 
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4-  M  ^  4- 
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i 


i 
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COMPUTER  CODING 
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SAMPLE  OUTPUT 


The  following  pages  give  the  program  output  for  the  example  case  presented  previously. 
This  output  is  fairly  self-explanatory,  although  the  following  points  are  worthy  of  men¬ 
tion. 

The  results  for  the  first  two  field  points  only  (11  were  requested)  are  included  for  the 
sake  of  brevity. 

It  should  be  remembered  that  the  program  calculates  the  sound  field  observed  when  the 
helicopter  rotor  is  positioned  at  the  nominal  position  X^,  Y^,  Z^.  This  sound  was 

actually  generated  by  the  rotor  when  it  was  in  some  other  position  (at  the  retarded 
time).  This  position  is  denoted  in  each  case  by  the  "retarded  hub  coordinates." 

Four  sound  harmonics  are  calculated  since  M  was  specified  as  36  so  that  (W/2&)  =  4. 5. 
The  Doppler  shift  effect  can  be  noted  in  the  slightly  different  frequencies  observed  at 
the  two  positions. 
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PROGRAM  .HERON  1 
CASE  H.1, 81, 

IN  LcVCL  TLliHT  at  40  KTS ,  OATa  FROM  NASA  TM.  X-952  TABLE  8, 


u6  JMcTR  Y  OF  ROTOR 

1  4 

BLADES 

2,4006  01 
2 , 1  /(It  02 

FT, 

R|P .M, 

radius 

1  ,  000b:  0(1 

FT. 

flap  hinge  offset 

1  ,  0 0 Ot  00 

FT. 

drag  hinge  offset 

0 

DEGREES 

azimuth  phase 

-1 , b’Oe-02 

DEGREES 

rearrard  shaft  inclination 

-1,8706-02 

DEGREES 

SHAFT  INCLINATION  to  starboard 

4, 7806  01 

FT/SEC 

VELOCITY  IN  X-DJRECTIQN 

0 

FT/SEC 

veloctiy  in  y. direction 

0 

FT/SEC 

velocity  IN  Z-DIHECTIOn 

6 , 7S0fc  01 

FT/SEC 

resultant  VELOCITY 

0 

FT. 

X  HUB  ORDINATE 

(1 

FT  . 

Y  HUB  ORDINATE 

2.0006  02 

FT, 

Z  HUB  ORDINATE 

dlAOc  lOaoINu  and  motion  OATa  (PARALLEL  AND  NORMAL  TO  THE  SHAFT 


NUMBER  OF  LOADING  STATIONS*  7 

NUMBER  OF  AZ  I  m  jth  INTEGRATION  POINTS*  34 

rORM  iJF  INPUT.  HARMONIC 


dLAUfe  FLAPPING  HARMONICS  (DEGREES) 


A  0 

3.307b  00 

A  1 

3 , 1 60b  -  0 1 

A  2 

1 , 400t-Q2 

A  i 

1,8006-01 

A  4 

7, 7006-02 

A  7 

A.nOOb-03 

A  6 

3.P006-02 

A  7 

1,4006-02 

a  a 

2,6006-02 

A  9 

6,0006*03 

A10 

1  ,  P  0  0  6  •  U  2 

li  1 

-4,7206-01 

a  2 

fl, 700E-02 

d  i 

*6,8006*02 

H  4 

8.8006-02 

b  8 

*4  ,  <5  00b*Q2 

d  6 

?.4QOb*02 

a  7 

ft, P00b-03 

b  8 

0 

H  V 

-4,000t-03 

HI  0 

?, C00E-03 

bLADt  LAGGING  HARMONICS  t  DEGREES ! 
A  0  -7,000b  00 
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".’MfcSSP'"  V**-' 


bLADE  ELEMENT  ELASTIC  DISPLACEMENTS  NORMAL  TO  THE  SHAFT 


RAOUL  station 

A  0 


7,000  11,200  15,400  21,000  23,800 

0  0  0  0  0 


25,200  26,600 

0  0 


bLADE  ELEMENT  ELASTIC  DISPLACEMENTS  PARALLEL  TO  THE  SHAFT 


RADIAL  STATION 
A  0 


11,200  15,400  21,000  23.800  25,200  26.600 

0  Of!  0  0  0 


blAllfc  THRUST  LOAD  HARMONICS  t  L  tip  > 


K*Di Al 
A  0 
A  1 
A  2 
A  3 
A  4 
A  5 
A  6 
A  7 
a  e 
A  9 
A10 
b  l 
b  2 
b  3 
b  4 
b  5 

b  6 

b  7 
b  d 
b  9 
blit 


STATION  7,000  11.200  15,400  21,000  23.600  25,200  26.600 

1 . 4  71 1  02  3.334E  02  5,66fifc  02  8,083fc  02  4.543E  02  3,0l4fc  02  3.081E  02 

-1.306E  01  3.024E  Of!  9.652E  01  1,3816  01-1. 779E  01-6, 166E  00  6.350E  00 

2.731E  01  5.690E  01  9,5566  01  1,2936  02  6,4496  01  2.340E  01-5.678E  00 

•b , 665t- 0 1  1.749E  01  3,1046  01  6,637t  01  7.7296  01  6.8766  01  6.693E  01 

5,9836  00-2,1176  00  7.877E  00  1,0566  01  8.8206-01  1.1426  01  3,7626  01 

6,4836  00  2,5206  00  5,8756-02-2,1826  01-1.8936  01-2,4016  01-2.9066  01 
-3.1576  00  4,9406  00  1,7636-01  1,7746  05  1,9406  00-1,2431:  00.2, 967E  01 
-5,5406  00-1,2046  01  1,2006  01  1,2706  01  4,3856  00  1,1046  01  1,0206  01 

-7,3136  00-3.3256  00-4,5286  00-1,9256  01  3.377F  00  2,0166-01  2,0906  01 

6,7606  00  7,8606  00-1, 570E  01*3,6286  00  3,7806-01-3.5626  00  1.428E  00 
.AilbOb  00-4,2326  00  1,0056  01  1,5676  01  5,4686  00  2,4706  00-1, 109E  01 

8,1156  05  9,5816  01  1,0856  02-3,7306  00-4,6126  01-4, 753E  01-4,0006  01 

1,9956  00  5,3406  00  1,0176  01  1,2256  01-2,0416  00-6,8546  00-1.403E  01 

2,1336  01  1,0576  00  1,3506  03-9,8786  00-2.2656  01-5,1586  00-1,3686  01 

6,593b  00-1,1796  01-1,2886  01-1,4116  01-2, 021E  01-1.4706  01  1.714E  00 

-2,5476  00-9,2236  00-1,5826  01-3,2866  01-1,0486  01-1,7911.  01-1,7396  01 

-5,6506  00-1,4266  U 1 -1 , 6826  01*1,3366  01-2.1046  Ul-7.3426  00-2.026E  01 

7. 7556-01-1, 7646  00-9,9356  00  6,5526  00  2,2936  0(1  3.226E  00  2,5206  0O 

3.7136  00  1,2636  00  4,5286  00-6,3006  00  1.021E  01  1,7646  00  1,4/26  01 

-'.,2166  on  1,7646  00-5, 232E  00  5,8466  00-1, 023E  01  1.1  /66  00  6,7876  00 

-1.4416  00-1,2636  00-9,2306  00  3 , 5286-01 -5 , 645E  00-1,7816  00-9,6106  00 


bl  A  1)6  DWAU  LOAD  HARMONICS  U.8S) 


HnDI AU  S  T  A  1  ION 

A  0 
A  1 
A  £ 


7,000 

0 

0 

II 

0 

0 

0 

0 

0 

ll 

0 

0 

0 

3 

0 

0 

0 

0 

0 

3 

0 

0 


11,200 

0 

U 

ll 

0 

u 

0 

0 

0 

u 

II 

u 

0 

0 

a 

u 

n 

0 

0 

3 

0 

0 


15,400 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


21,000 

0 

n 

n 

n 

n 

n 

0 

0 

0 

0 

0 

0 

n 

0 

0 

n 

n 

0 

0 

0 

0 


23,8011 

0 

0 

0 

0 

n 

i) 

0 

0 

0 

0 

0 

0 

0 

0 

c 

0 

0 

0 

0 

0 

0 


25,200 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

c 

0 

0 

0 

0 

0 

0 


26,600 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 
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COMPUTED  SOUND  FIELD 
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Th»e  report  contains  a  comprehensive  description  of  a  computer  program  developed  for  th>» 
numerical  evaluation  of  the  helicopter  noise  equations  derived  in  USAAVLABS  TR  68-60 ,  "Studies 
of  Helicopter  Rotor  Noise."  It  is  completely  self-contained  in  that  the  program  details  ore 
described,  starting  from  Iwo  basic  acoustic  equations  and  covering  methods  by  which  these 
equations  are  applied  to  I  he  rotor  noise  problem.  Program  flow  diagrams  and  a  complete  listing 
are  presented  together  with  input  instructions  and  sample  inputs  and  outputs.  The  program  is 
written  in  FORTRAN  IV for  the  CDC  3300  Computer. 
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